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This  report  presents  an  analysis  and  experimental  data  in  its 
support  for  a  laminar  jet  angular  rate  sensor  (LJARS) .  An  LJARS 
is  a  device  in  which  a  Coriolis  force  causes  a  jet  to  deflect  as 
a  function  of  applied  angular  rate.  This  jet  deflection  is  sensed 
as  a  differential  pressure  or  flow.  Both  the  theory  and  the  data 
show  that  the  differential  pressure  is  a  linear  function  of  applied 
angular  rate.  The  report  considers  the  fluid  mechanic  and  fluid  — , 
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dynamic  aspects  of  the  device  and  discusses  in  detail  such  areas 
as  jet  pressure  recovery,  dynamic  response,  and  transport  delay. 
In  addition,  the  report  expands  considerably  on  the  theory  by 
presenting  numerous  empirical  guidelines  for  the  design  and  the 
implementation  of  these  devices.  Included  is  a  brief  heuristic 
argument  for  the  determination  of  the  transition- to- turbulence 
point  and  methods  for  effecting  suitable  temperature  and  null 
offset  compensation.  Environmental  temperature  and  pressure 
effects  are  detailed.  Finally,  typical  LJARS  characteristics 
are  summarized  and  presented.  The  important  characteristic 
parameters  are  thresholds  of  less  than  0.1  deg/s,  offsets  of 
less  than  0.5  deg/s,  and  sensitivities  of  0.04  Pa/ (deg/s).  The 
LJARS  output  is  easily  amplified  by  laminar  proportional  ampli¬ 
fiers.  Other  desirable  characteristics  of  the  LJARS  are  an 
estimated  cost  of  less  than  $100,  high  reliability,  and  instant- 
on  performance. 
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FOREWORD 


Since  their  inception  in  the  early  1970's*  the  development  of 
laminar  fluidic  control  components  has  been  documented  almost  solely 
with  regard  to  the  laminar  proportional  amplifier,  yet  the  laminar  jet 
angular  rate  sensor  has  been  as  valuable  as  a  control  component,  if  not 
more  valuable.  This  report  attempts  to  remedy  this  imbalance  of  techni¬ 
cal  information  by  compiling  most  of  the  work  that  has  been  found  in 
laboratory  notebooks  of  various  researchers.  Some  data  were  compiled  at 
the  Harry  Diamond  Laboratories  (HDL)  by  Hebert  Foster  (formerly  of  HDL), 
Steve  Tenney,  Charles  Paras,  George  Mon,  and  Kenji  Today  and  Margaret 
Hackert  and  Kevin  Hippie,  former  summer  employees  of  HDL.  Other  data 
have  been  obtained  by  HDL  on  contract  with  industrial  concerns : 
AiResearch  Manufacturing  Co .  of  Arizona;  McDonnell  Douglas,  and  TriTec, 
Inc.  The  analysis  for  the  most  part  is  of  the  authors. 

It  is  hoped  that  this  report  will  lay  the  technical  foundation  for 
systematic  sensor  design. 
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1 .  INTRODUCTION 


Fluidic  angular  rate  sensing  is  older  than  fluidics  itself.  Born 
out  of  the  need  to  replace  costly  and  unreliable  mechnical  rate  gyro¬ 
scopes,  fluidic  rate  sensors  have  carved  themselves  out  a  niche  in  this 
area.  Perhaps  the  most  familar  fluidic  sensor  is  the  vortex  rate  sensor 
(VRS).  This  device  makes  use  of  the  induced  weak  vortex  on  an  inward 
radial  flow  to  provide  a  signal  as  a  function  of  applied  rate. 

Numerous  technical  papers  discuss  this  device,  but  perhaps  most 

relevant  and  useful  are  those  by  Ostdiek,1  Kirshner,2  and  Hedeen.3  The 

first  two  discuss  the  analytical  aspects  in  detail ,  and  the  last  de¬ 
scribes  the  use  of  such  sensors  in  systems.  More  recently,  Ringwall*  has 
described  a  device  called  the  Fluidise,  which  uses  a  flow  of  liquid  in  a 
helical  set  of  passages  to  replace  the  spinning  mass  of  a  gyroscope. 
Rate  is  detected  by  a  force  balancing  system  on  the  helical  flow  con¬ 
tainer.  This  device  has  a  moving  part.  Its  performance  is  discussed  by 
Chen  and  Neradka4  in  a  report  comparing  various  fluidic  rate  sensors. 

The  concept  of  the  laminar  jet  angular  rate  sensor  (LJARS)t  was 

proposed6  in  the  early  1940's,  but  the  patents  for  jet  angular  rate 

sensors  appeared  much  later  (about  1970)  and  are  held  by  Reader6  and 


1A.  Ostdiek,  viscous  Vortex  Rate  Sensor,  Harry  Diamond  Laboratories 
HDL-TR-1555  (November  1971). 

2J.  M.  Kirshner,  A  Survey  of  Sensors — Part  II,  Proceedings  of  HDL 
Fluidic  State-of-the-Art  Symposium,  II,  Harry  Diamond  Laboratories 
(October  1974). 

3J.  O.  Hedeen,  Performance  and  Application  of  the  Vortex  Rate  Sensor, 
Proceedings  of  HDL  Fluidic  State-of-the-Art  Symposium,  II,  Harry  Diamond 
Laboratories  (October  1974). 

4K.  Chen  and  V.  Neradka,  Development  of  Critical  Components,  TriTec, 
l.JC.,  Columbia,  MD,  DAAG-39-77-C-0173  (July  1978). 

5ff.  Ziebolz,  Characteristics  of  Hydraulic  and  Pneumatic  Relays  as 
Energy-Converting  Devices,  Instruments,  _15_  (September  1942). 

6T.  Reader,  Laminar  Fluidic  Devices,  U.S.  Patent  3,662,772 

(16  May  1972). 

*C.  G.  Ringwall,  Summary  Report  on  Single  Axis  Stabilization  System, 
General  Electric  Co.,  Schenectady,  NY,  DAAG-39-76-C-0137  (February 
1977). 

tThis  device  is  also  known  as  LARS,  for  laminar  angular  rate  sensor,- 
however,  since  most  high-performance  vortex  rate  sensors  also  are 
laminar  and,  hence,  laminar  angular  rate  sensors,  "LJARS”  is  used  to 
avoid  confusion. 
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Moore.7  The  way  that  this  device  works  is  that,  as  a  jet  leaves  a 
nozzle,  it  tends  to  travel  in  a  straight  line  in  inertial  space — yet  the 
device  structure  is  rotating  as  shown  in  figure  1.  By  the  time  that  the 
particles  get  across  the  device,  the  structure  has  moved,  and  the  jet 
deflection  is  relative  to  the  structure.  This  relativity  causes  a 
pressure  difference  on  outputs  as  in  an  LJARS  or  a  temperature  differ¬ 
ence  on  thermistors  as  in  a  "Super  Jet."8  The  pressure  or  temperature 
difference  is  related  directly  to  angular  rate. 


TIME  HISTORY  OF  ONE 
PARTICLE  AND 
SPLITTER  OVER 
0  1  t  S  3 

DIRECTION  OF 
PARTICLE  LEAVING 
ORIGIN  POINT 
AT  t  =  3  AND  LOCUS 
OF  INTERMEDIATE 
PARTICLES 


POSITION  OF 
SPLITTER 
AND  PARTICLE 
t  =  0 


POSITION  OF 
SPLITTER  AND 
PARTICLE  AT 
t  =  3  ; 


Figure  1.  Coriolis  jet  deflection. 


7A.  G.  Moore,  Angular  Movement  Sensing  Device,  U.S.  Patent  3,500,691 
(17  March  1970). 

8T.  M.  Rankin,  A.  G.  Moore,  and  M.  C.  Schuemann,  Fluidic  Angular  Rate 
Sensor — A  Replacement  for  Rate  Gyroscopes?  Johns  Hopkins  Applied  Physics 
Laboratory  Technical  Digest  (March-April  1969). 
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Research  on  the  LJARS  was  spurred  by  the  limitations  on  the  VRS, 
namely,  the  high  power  (or  flow)  consumption,  the  very  high  output 
impedance,  and  the  relatively  low  bandwidth.  Ihe  high  power  require¬ 
ment4 * * * * * 10  of  the  VRS  tended  to  place  a  burden  on  existing  power  supplies, 
especially  if  hydraulic  oil  was  used.  The  high  output  impedance  effec¬ 
tively  reduced  the  high  sensitivity  in  circuits  where  the  signal  was 
matched  to  other  fluidic  components.  Most  devices  exhibited  output 
impedances  one  to  three  orders  of  magnitude  greater  than  the  input  of  a 
fluidic  circuit;  hence,  the  net  sensitivity  was  decreased  by  the  same 
order  of  magnitude.  The  bandwidth  is  limited  by  the  time  required  to 

fill  the  relatively  large  volume  of  the  device. 

In  contrast,  an  LJARS  requires  very  little  power  to  establish  a  rate 
sensitive  jet,  has  the  output  impedance  of  a  fluidic  amplifier  (hence, 
should  match  well)  ,  and  has  a  bandwidth  limited  only  by  the  fluid 
transport  time  across  a  small  gap.  Early  attempts  at  analyzing  and 
using  the  LJARS9-11  led  to  devices  shown  in  figure  2.  A  long  nozzle 
was  used  because  it  was  believed,  at  that  time,  that  this  length  would 
promote  a  greater  laminar  range. 10,11  Subsequently,  it  was  found  that 
short  nozzles  appropriately  designed  would  provide  greater  pressure 
recovery  and  gain  in  laminar  amplifiers  (hence,  greater  LJARS  sensitiv¬ 
ity)  and  maintain  superior  laminar  range.12'13  The  early  devices, 
however,  did  prove  an  important  point  in  that  they  demonstrated  the 
feasibility  of  laminar  jet  angular  rate  sensing  and  its  advantages, 

including  a  remarkable  (for  that  time)  dynamic  range  in  excess  of 

14,000.  Dynamic  range  is  defined  as  the  maximum  input  signal  (that 
required  to  saturate)  divided  by  the  threshold  or  smallest  input. 
Young11  reports  an  LJARS  with  a  threshold  of  0.05  deg/s  and  a  maximum 
rate  of  700  deg/s — hence,  a  dynamic  range  of  ~14,000. 

4AT.  Chen  and  V.  Neradka,  Development  of  Critical  Components,  TriTec , 

Inc.,  Columbia,  MD,  DAAG-39-77-C-0173  (July  1978). 

9T.  Reader,  T.  Shaffer,  J.  Goren,  and  R.  Kantola,  Feasibility 

Investigation  of  a  Laminar  Rate  Sensor  (LARS),  Phase  I  Report,  General 

Electric  Co.  Reentry  and  Environmental  Systems,  Schenectady,  NY,  69SD698 

(June  1969). 

10P.  L.  Jacobs,  Theoretical  Analysis  of  a  Two-Dimensional  Laminar  Jet 
Rate  Sensor,  Proceedings  of  HDL  Fluidic  State-of-the-Art  Symposium,  II, 
Harry  Diamond  laboratories  (October  1974). 

Hr.  Young,  Development  of  a  Laminar  Angular  Rate  Sensor  (LARS),  Pro¬ 
ceedings  of  HDL  Fluidic  State-of-the-Art  Symposium,  II,  Harry  Diamond 
Laboratories  (October  1974). 

12f.  Manion  and  T.  Drzewiecki,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I_,  Harry  Diamond  Laboratories  (October  1974). 

1 3T.  M.  Drzewiecki,  D.  N.  Wormley,  and  F.  M.  Manion,  Computer-Aided 
Design  Procedure  for  Laminar  Fluidic  Systems,  J.  Dyn.  Syst.  Meas.  Con¬ 
trol,  97,  Series  G  (December  1975). 
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Figure  2.  Early  laminar  jet  angular 
rate  sensor  design. 


Rankin  et  al8  describe  the  axisymmetric  Super  Jet  WARS  as  a  viable 
replacement  for  rate  gyroscopes,  quoting  such  performance  parameters  as 
drift  of  0.001  deg/s  per  day  and  resolution  of  less  than  0.002  deg/s. 
These  are  remarkable  even  for  good  rate  gyroscopes,  considering  that 
earth  rate  at  the  equator  is  0.004  deg/s. 

Hiis  report  lays  the  analytical  groundwork  for  the  coherent  design 
of  a  planar  LJARS.  The  report  is  divided  into  two  sections,  (1)  the 
purely  theoretical  analysis  verified  by  experimental  data  and  (2)  the 
design  and  applications  guide.  This  second  section  is  subdivided  into 
critical  areas  that  relate  to  the  generation  of  a  real-life  system  to 
include  operating  range,  temperature  (environmental)  sensitivity,  shock 
and  vibration  sensitivity,  null  offset,  drift,  dynamic  range,  power 
supply  considerations,  noise  effects,  operating  fluids,  fabrication,  and 
testing.  Where  feasible,  analytic  tools  are  used  to  predict  perform¬ 
ance.'  *  Where  ‘not, '  rules  of  thumb  and  the  experience  of  the  authors  and 
others  are  presented. 


2.  BASIC  ANALYSIS 

2. 1  Free  Jet  Deflection  in  Rotating  Field 

The  basic  phenomenon  of  the  WARS  relies  on  the  time  that  it 
takes  a  particle  leaving  the  nozzle  to  reach  the  sensing  region,  as 
illustrated  in  figure  1.  During  that  time,  the  outputs  or  sensors  have 
moved  relative  to  their  initial  locations.  As  a  result,  the  deflection 
of  the  jet  is  relative  to  the  outputs  (jet  deflection  due  to  Coriolis 
forces).  With  a  rotating  enclosed  body  like  a  fluidic  amplifier,  the 
field  into  which  the  jet  enters  is  rotating  with  the  structure  due  to 
the  shear  exerted  on  the  fluid  by  the  solid  boundary  planes.  This 
rotation  reduces  the  total  jet  deflection  over  that  which  would  have 


8T.  M .  Rankin,  A.  G.  Moore,  and  W .  C.  Schuemann,  Fluidic  Anguleur  Rate 
Sensor — A  Replacement  for  Rate  Gyroscopes?  Johns  Hopkins  Applied  Physics 
Laboratory  Technical  Digest  ( March-April  1969). 


1 


rir 


L2 


occurred  if  the  field  had  been  stationary.  Consider,  therefore,  the  net 
displacement  of  a  fluid  particle  traveling  with  a  velocity  ucg  in  a 
radial  direction  (relative  to  the  axis  of  rotation),  with  a  superimposed 
solid  body  rotational  velocity,  as  shown  in  figure  3. 


SPLITTER 


Figure  3.  Geometry  and  coordinates  for  rotational 
deflection. 


The  particle  displacement,  y. ,  at  a  distance,  x,  from  the 
origin  is  the  net  tangential  motion  at  tangential  (lateral)  velocity,  v, 
over  the  period  of  time,  t,  that  it  took  the  particle  to  reach  x  from 
the  origin.  At  the  average  jet  velocity,  the  lateral  velocity  is 

v  =  r9  cos  Y  , 

where  the  angle  Y  =  tan-1  yj/r  so  that,  for  transit  velocities  much 
greater  than  the  lateral  velocity,  cos  y  a  1.  Therefore,  the  elemental 
lateral  displacement  is 


dy.  =  r9  dr 
D 
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■ito  -  f ’-*■  iMim  l  l  V  -  ;• 


(1) 


and  the  elemental  transit  time,  di,  Is 


dt  = 


dx 


(2) 


so  that  the  total  lateral  displacement  of  the  particle  in  the  field  is 


fr'  Fsp  .  ax 

ri-Jo  y'X  r6^‘ 


(3) 


But  r  =  rn  +  x- so  that 


[*8P  r  +  x 

’  .  =  9  I  — -  dx  , 

3  Jo  »cZ 


(4) 


where  x  is  the  distance  from  the  nozzle  to  the  splitter  as  in  fig¬ 
ure  3.  fee  distance  traversed  by  the.  splitter  (or  the  outputs)  at  the 
tangential  velocity,  rsp0  *  (rn  +  xsp)0  i**  the  total  transit  time. 


=  [*SP  <*. 
Jo  V 


is 


y  =  /r  +  x  \  9 
xsp  In  spl 


Jo  UcZ 


(5) 


The  displacement,  yr,  of  the  jet  relative  to  the  structure  is  therefore 
the  difference  of  these  two  displacements, 


y  =  y  -  y.  = 
r  xsp 


,.(,♦,)»  f*sp — -  e  h 

3  \  "  Bpf  J0  Mcl  Jo  ncl 


dx  • 
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The  terms  containing  the  center  of  rotation,  rn,  drop  out  and,  hence 4 


.  r 

'r  ‘°Jo 


sp  X  -  X 

S2 -  . 

uc* 


(6) 


For  the  high  Reynolds  number  case  where  the  center-line  velocity  is 
constant,  equation  (6)  reduces  to  the  familiar  form 


0x2 


J!£ 


2u 


cl 


(7) 


where  uQ ^  =  Qs/ob^,  where  Qs  is  the  jet  flow  rate,  bg  is  the 

nozzle  width,  and  o  is  the  aspect  ratio,  h/bg. 

It  now  remains  to  identify  the  cases  in  which  the  average 
center-line  velocity,  ucj,,  is  not  constant.  At  low  Reynolds  numbers, 
the  two-dimensional  laminar  jet  decays  as  the  one-third  power  of  the 
distance  from  the  virtual  origin  so  that 


u 

cl 


(8) 


The  virtual  origin,  x^,  is  determined  by  matching  the  jet 
center-line  velocity  to  that  of  the  exit  of  a  nozzle  with  fully 
developed  (parabolic- like)  flow.  The  reasoning  behind  this  choice  is 
that  the  parabolic  profile  and  the  hyperbolic  secant  profile  are  similar 
near  the  center,  as  is  shown  by  Jacobs.10  The  actual  center-line 
velocity  is 


u 


cl 


0.4543 


(9) 


10P.  L.  Jacobs,  Theoretical  Analysis  of  a  Two-Dimensional  Laminar  jet 
Rate  Sensor,  Proceedings  of  HDL  Fluidic  State-of-the-Art  Symposium,  II , 
Harry  Diamond  Laboratories  (October  1974).  ' 


where 


and 


Uj 


1/2 


K  = 


2cflb  P 
0  S  S 


where  K  is  the  kinematic  raomentun  flux,  v  is  the  kinematic  viscosity, 
is  the  volumetric  flow  rate  of  the  supply,  Pg  is  the  supply  pressure, 
is  the  density,  and  Cg  is  the  momentun  flux  discharge  coefficient. 
Therefore, 


or 


where 


vo 


0. 4543  I3 


k£ 

v 


0.01172  —  N 

3  R 
r  * 


(10) 


c9  ■  ’-,5i 


(based  on  a  parabolic  fully  developed  profile  from  a  rectangular 
nozzle).  Therefore, 


X  =  0. 01 55c .N_  . 
vo  d  R 


(11) 
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to  a 


The  center-line  velocity  is  assumed  to  decay  as  in  equation  (12): 


,(?) 


(.  ♦ 

li  +  =*-Y 

\  xvo/ 

\  VO  / 

1/2 

\l/3 


It  has  become  obvious  that  three  different  functional 
relationships  govern  the  average  jet  velocity,  in  three  Reynolds 

number  regimes  where  the  potential  core  (the  distance  over  which  the 
center-line  velocity  is  constant)  determines  the  velocity. 


jfff  ■ 


high  NR,  xCQre  2  xflp, 


JW 


1/2 

\l/3 


x’  "  x  **  xcore'  0  *  xcore  <  xsp' 


S' 


,  low  Nr,  xcore  ^  0. 


Hie  core  length  concept  is  illustrated  in  figure  4.  The  core 
length  is  defined  by  the  point  where  the  nozzle  boundary  layers  merge, 
or  5/bs  =  0.5.  Using  this  criterion  and  the  fact  that  the  boundary 
layer  grows  as  the  square  root  of  downstream  distance,  we  obtain 


(Xn  +  xcore] 


n 


core 


m 


and 


1  il 

1 5 

\2 

4  b 

-) 

L  \ 

s)  exit 

where  3^  is  the  effective  length  of  flat  plate  required  to  generate  a 
nozzle  exit  boundary  layer  thickness  of  6/balex^t. 


BOUNDARY  LAYER 


Figure  4.  Physical  interpretation  of  virtual  origin 
and  potential  core. 


Hie  boundary  layer  thickness  is  related  to  the  displacement 
thickness,  6*,  which  is  known  from  the  discharge  coefficient, 12 


£  -  “-P  -  -  (’  -  «a)  S]1'2 


(15) 


12F.  Manion  and  T.  Drzewiecki ,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  _I,  Harry  Diamond  Laboratories  ( October  1974). 


For  the  zero  pressure  gradient14'15  at  the  nozzle  exit,  6/6*  *  10/3  and 
6/bs  *  5(xn/NR)1'2;  therefore,  the  effective  nozzle  length  becomes 


X  = 
n 


(16) 


From  observations  of  flow  visualization  experiments,  it  is 
found  that  the  core  length  is  considerably  longer  than  that  given  by 
equation  (14),  probably  because  the  bounded  jet  configuration  tends  to 
reduce  entrainment  effects  and,  hence,  does  not  degrade  the  center 
velocity  as  much.  When  one  choses  the  fictitious  conditions  for 
boundary  layer  thickness  at  the  end  of  the  potential  core  as  VhaIcore  = 
1.19  instead  of  0.5  as  originally  assumed,  the  resulting  equation  gives 
good  results  in  the  overall  algorithm.  Hence,  equation  (14)  becomes 


r  i 


The  general  form  for  the  jet  deflection  thus  is 


b 


s 


-  x)  f(X)  dx  , 


(18) 


where  Xgp  is  now  the  splitter  distance  normalized  by  the  nozzle  width 
and 


14T.  M.  Drzewiecki ,  Interpretation  of  Surface  Static  Pressure  Distri¬ 
butions  in  Fluid  Amplifier  Applications ,  Harry  Diamond  Laboratories 
HDL-TR-1627  (July  1973). 

15H.  Scblichtirtg ,  Boundary  Layer  Theory,  McGraw-Hill  Book  Co.,  New 
York  (1960). 


19 


general  case  when  0  <_  xcore  <  x  degenerates  into  the  two 
cases  when  at  the  limits;  hence, 


other 


and,  performing  the  nontrivial  integration. 


20 


V 


X _ _  >X„  ) 

\  core—  spj 


2c . 


(3 


■  X2 

1/2  8P 


(20) 


which  is  identical  with  equation  (7) . 
Equation  (19)  reduces  to 


*r 

b_  " 


bs6 


cd\  p 


>  V  1 
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)a_  X2 

j  28  xvo 

-  4  XvoXsp| 

( 

) 

(21) 


for  Xcore  =  °»  that  is,  when  there  is  no  core  length.  The  relationship 
for  deflection  can  then  be  written  as 


(22) 


where  ZO  represents  one-half  the  square  of  the  effective  distance  from 
nozzle  to  splitter.  When  the  core  length  xcore  >_  x  ,  ZO  =  x|  /2.  The 
coefficient  of  ZO  may  be  rewritten  as 

V 

Cd\  P  / 

Now  the  output  pressure  is  related  to  deflection12  by 


12F.  Manion  and  T.  Drzewiecki,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I_,  Harry  Diamond  Laboratories  (October  1974). 
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(24) 


where  Cj^gg  is  the  percentage  of  the  momentum  lost,  BQ  is  the  normalized 
width  of  the  output  channel,  and  Prec  is  the  jet  centered  recovered 
pressure  at  blocked  conditions. 


We  define  normalized  sensitivity  as 


(25) 


To  get  APq/0  in  pressure  units  per  degree  per  second,  we  multiply  S  by 
Pgb^/57. 3v  so  that 


AP 
_ o 

• 
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_2 _ 

57.3 


P 

rec 


P 

s 


zo  . 


(26) 


Sensitivity  becomes  dimensionless  when  divided  by  absolute  viscosity, 
p.  This  quantity,  APQ/p0  ,  is  not  the  normalized  sensitivity,  but  is  a 
dimensionless  sensitivity. 


__2 _ 

57.3 


P 

rec 


P 

s 


ZO  . 


(27) 


To  obtain  sensitivity,  one  need  only  to  multiply  the  dimensionless 
sensitivity  by  the  viscosity  of  the  fluid  being  used.  The  customary 
units  are  pei/( deg/s);  SI  units  are  Pa/( deg/s).  The  sensitivity  of  a 
device  operating  in  a  given  median  is  independent  of  actual  size  and 
depends  only  on  the  absolute  viscosity  of  the  fluid.  The  difference  in 
operation  between  LJARS's  of  different  sizes  is  in  maximan  rate  since  a 
larger  device  saturates  before  a  small  one  due  to  the  longer  transit 
time. 
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2.2  Augmented  Jet  Deflection  in  Presence  of  Controls 


When  a  jet  deflection  not  due  to  a  control  signal  occurs  in  a 
proportional  amplifier  configuration,  a  pressure  difference  develops 
across  the  jet  and  affects  the  jet.  Apparent  or  real  deflections  occur 
due  to  misalignments,  rate  inputs,  and  geometric  anomalies.  In  a  rate 
sensor,  therefore,  the  deflection  can  be  augmented  by  this  positive 
feedback,  which  is  illustrated  in  figure  5. 


Figure  5.  Deflection  augmentation  due  to  offset. 


The  physical  explanation  for  this  phenomenon  is  straightfor¬ 
ward.  Since  the  jet  is  closer  to  one  control  edge  them  to  the  other, 
the  resistance  of  that  clearance  is  higher.  So  more  entrainment  flow 
must  be  supplied  fran  the  control.  Since  control  pressure  is  constant 
and  equal  on  both  controls,  more  control  flow  requires  that  the  jet  edge 
pressure  be  lower.  Conversely,  the  reduced  resistance  of  the  wider  side 
increases  that  pressure.  This  net  pressure  deflects  the  jet  additional¬ 
ly.  Deflection  stops  when  the  centripetal  forces  of  the  curving  jet 
balance  the  applied  pressure  difference. 

To  compute  this  additional  deflection,  we  may  use  the  block 
diagram  systems  approach.  An  input  deflection  due  to  rate  at  the  end  of 
the  control  region  causes  a  control  pressure  difference,  which  causes 
more  deflection,  which  increases  the  deflection  at  the  splitter 
(fig.  6). 


LATERAL 

DEFLECTION 


ANGULAR  DISTANCE  PRESSURE 


Figure  6.  Positive  feedback  for  added  deflection  due  to  presence  of 
controls. 


Section  2.1  has  shown  how  the  transfer  function  is  deter¬ 
mined.  This  section  derives  the  transfer  functions  for  induced  pressure 
differential  and  for  induced  deflection. 

2.2.1  Pressure  Caused  by  Offset 

In  figure  7,  consider  the  control  volume  near  the  control.12 

The  application  of  the  incompressible  continuity  equation 

results  in 


12F.  Manion  and  T.  Drzewiecki,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I_,  Harry  Diamond  Laboratories  (October  1974). 
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=  0  , 


(28) 


I 


where 


2  .  +  8  .  +  Q 
cl  vl  y 


enet 


Qv1 

Qy 

®enet 


control  flow, 

flow  from  vent  through  undeflected  jet-edge  space 

that  has  resistance  Ry, 

flow  through  deflected  space, 

net  entrained  flow,  that  is,  entrained  flow,  Qe, 
minus  spill  flow,  QSB.  (Part  of  entrainment 
requirement  is  self-satisfied  by  spill  flow.) 


The  control  resistance,  Rc,  is  known,  as  is  the  resistance, 
Ry,  of  the  centered  jet, 12  where  resistance  is  defined  R  =  AP/Q. 
Equation  (28)  can  thus  be  written  in  terms  of  the  control  volume  pres¬ 
sure. 


P 


cl 


Pj1 


+ 


(Ce  "  Csb)  =  0  '  <29> 


PRESSURE  LATERAL 


12F.  Manion  and  T.  Drzewiecki ,  Analytical  Design  of  laminar  propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I_,  Harry  Diamond  Laboratories  (October  1974). 

16r.  M.  Drzewiecki ,  Fluerics  38:  A  Computer  Aided  Design  Analysis  for 
the  Static  and  Dynamic  Port  Characteristics  of  Laminar  Proportional 
Amplifiers,  Harry  Diamond  Laboratories  HDL-TR-1758  (June  1976). 


where  the  flow  through  the  deflected  space  is  treated  as  an  orifice  flow 
through  an  area  yBCh  with  unity  discharge  coefficient.  When  equation 
(29)  is  normalized  by  the  supply  jet  flow  Qs  =  c(Jbsh(2Ps/p)1' 2  ,  it 
becomes  as  shown  in  equation  (30),  in  which  the  prime  (')  notation 
indicates  quantities  normalized  by  the  corresponding  supply  parameter. 
For  example,  P‘  =  P/Ps,  =  Qe/Qg,  «nd  Rt  =  H/V  where  Rs  =  Ps/Qs* 


Pc1 


P  1 

Pj1 


Rc 


P»  -  p  * 

V  pji 

+  - -  + 


^BC 


(K  -  "),) 


l/ 2  _ 


q;  +  =  0 


(30) 


Wte  now  drop  the  prime  notation  for  simplicity.  To  determine 
the  change  of  the  jet  edge  pressure  with  a  change  in  deflection,  we 
require  the  derivative  of  Pj1  with  respect  to  yr.  Ihus,  we  differ¬ 
entiate  equation  (30)  with  respect  to  yr  and  evaluate  at  yr  =  0  (the  jet 


centered  condition).  For  a  rate  sensor,  the  control  pressure 
remains  constant.  For  all  cases,  normally,  Pv  also  is  cc 
entrainment  flow  also  is  assumed  to  remain  constant  so  that 
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RcRv 
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(31) 


Similarly  for  side  2,  the  same  equation  is  written 
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Now  that  an  expression  has  been  obtained  for  the  change  of  the  jet  edge 
pressure,  it  remains  only  to  find  the  differential  applied  pressure  on 
the  jet.  The  pressure  at  the  edge  of  the  jet  can  be  found  by  expanding 
the  pressure  in  a  Taylor  series  in  deflection  and,  for  small  deflec- 


first 

two  terms .  Hence 

(yr)  ' 

9Pj1 

=  Pj1(0)  +  yr  3y 
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(-yr) 
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Since  Pj.,(0)  *  Pj2<°>* 
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or  the  transfer  function  for  jet  edge  differential  pressure  is 
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(36) 


itanion  and  Drzewiecki ^2  define  twice  the  change  in  spillback 
flow  with  respect  to  deflection  as  the  net  entrainment  coefficient,  a  ^ , 
which  is  calculated  from  the  parameters  of  the  problem.  An  improved 
expression  is  given  by  Drzewiecki.***  Therefore, 


3Q 


SB 


ay 


y=0 


(37) 


If  we  consider  Pv  =  0  so  that  all  measurements  are  relative  to  vent 
pressure,  Pv,  it  then  remains  only  to  find  an  expression  for  P  j(0)  . 
Manion  and  Drzewiecki12  do  this  by  applying  continuity  to  the  control 
leg  of  the  device  as  in  equation  (30)  at  y  =  0  so  that 


P.(0) 

3 


(38) 


When  these  results  are  substituted  into  equation  (36),  the  dc  transfer 
function  for  jet  pressure  is 


AP. 


y 


(39) 


12f.  Manion  and  T.  Drzewiecki ,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers ,  Proceedings  of  HDL  Fluidic  state-of-the-Ait 
Symposium,  Harry  Diamond  Laboratories  (October  1974). 

16T.  M.  Drzewiecki,  Fluerics  38:  A  Computer  Aided  Design  Analysis  for 
the  Static  and  Dynamic  Port  Characteristics  of  Laminar  Proportional 
Amplifiers,  Harry  Diamond  Laboratories  HDL-TR-1758  (June  1976). 
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One  observes  now  from  the  derivation  that  the  quantity  in  brackets  in 
equation  (39)  must  always  be  positive  so  that  when  the  control  flow, 
Pc/J^j,  exceeds  the  net  entrainment  flow,  the  effect  is  a  retarding  one, 
as  is  that  of  the  spill  flow.  Hence,  for  a  rate  sensor  with  close-in 
controls,  the  effects  of  the  controls  can  be  eliminated  by  judicious 
application  of  biasing  control  flow  so  that  AP-/y  =  0.  This  is  an 
important  point  and  is  discussed  further  in  section  3.4  on  null  offset 
compensation.  Equation  (39)  is  rewritten  to  accommodate  an  always 
positive  radical,  shown  in  parentheses  in  equation  (40). 


(40) 


It  is  clear  that  there  is  a  value  of  Pc  for  which  AP  /y  =  0  and  that, 
when  that  value  is  exceeded,  the  effect  is  a  net  reduction  in  deflec¬ 
tion. 


The  transfer  function  that  determines  the  pressure  caused  by 
an  offset  is  given  by  equation  (40). 

2.2.2  Augumented  Deflection 

As  a  result  of  the  differential  pressure  imposed  on  the  jet, 
the  jet  will  deflect  in  a  circular  arc  within  the  control  region  and 
have  a  deflection  at  the  end  of  the  controls  as  described  by  ftenion  and 
Mon*7  and  Manion  and  Drzewiecki*2  and  given  in  equation  (41). 


AP 


c 

4C. 


(41) 


where  Bc  =  b^bg  (the  normalized  control  width) . 


12F.  Manion  and  T.  Drzewiecki ,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I_,  Harry  Diamond  Laboratories  (October  1974), 

17F.  M.  Manion  and  G,  Mon,  Fluerics  33:  Design  and  Staging  of  Laminar 
Proportioned  Amplifiers,  Harry  Diamond  Laboratories  HDL-TR-1608 
(September  1972), 
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The  net  added  deflection  is  determined  by  the  transfer 
function  of  the  positive  feedback  loop  that  has  unity  feedback.  The 
gain  of  such  a  loop,  G^,  is 


G  =  __g _ L_ 

GL  1  -  G  _1  _  1 
G 


where  G  is  the  forward  gain;  hence, 


yBC  &Pj 


y__  AP. 
BC  2 


When  equations  (40)  and  (41)  are  substituted  into  equation 
(42)  and  the  result  is  rearranged,  the  result  is 


4Cn I R_  + 


enet  Rc 


The  jet  would  flow  in  a  straight  line  to  the  splitter  due 
only  to  control  deflection;  hence,  the  additional  displacement  at  the 
splitter  is  that  linear  contribution  from  equation  (43)  multiplied  by 
the  geometric  gain  factor  given  by  Manion  and  Drzewiecki.12 


-  1  . 


12 F. Manion  and  T.  Drzewiecki,  Analytic^  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I,  Harry  Diamond  Laboratories  (October  1974). 


The  last  transfer  function  required  is  the  deflection  at  the 
control  edges  due  to  rate.  For  this,  we  use  equation  (20)  with  the 
control  width  substituted  for  the  distance.  This  is  reasonably  accurate 
since  right  near  the  nozzle  exit  one  never  expects  there  to  be  much 
velocity  decay.  Thus, 


The  transfer  function  for  the  total  dimensional  deflection  is 
therefore  obtained  front  equations  (22)  and  (43)  to  (45)  as  is  expected 
from  the  block  diagram  of  figure  6. 


eq  (45)  (43)  (44) 


4c„  R  + 


2.3  Loss  of  Momentum  Due  to  Viscous  Friction 


The  loss  of  momentum  due  to  friction  on  the  bounding  planes  has 
been  addressed  by  Manion  and  Drzewiecki . 12  The  analysis,  however,  over¬ 
estimated  these  losses  at  low  Reynolds  numbers.  The  purpose  of  this 
section  is  therefore  to  remedy  this  loss  in  the  light  of  possible 
operation  of  laminar  devices  in  the  low  NR  region.  This  loss  occurs 
most  often  with  devices  with  very  long  nozzles  as  it  occurred  with  the 
early  LJARS.  Operation  with  a  long  nozzle  is  usually  undesirable 
because  of  the  dramatically  increased  sensitivity  to  temperature 
changes.  Nevertheless,  it  is  of  academic  interest  to  be  able  to  predict 
such  performance,  and  it  is  important  to  be  able  to  predict  an  LJARS 
performance  over  the  entire  laminar  range.  With  amplifiers,  one  is  not 
particularly  interested  in  the  region  where  the  gain  is  less  than  1,  but 
in  the  LJARS  there  is  always  rate  sensitivity  so  that  all  Reynolds 
numbers  are  of  interest. 

The  net  change  of  momentum  can  be  estimated  if  the  velocity 
distribution  at  two  stations  is  known.  The  momenta  when  the  jet  exits 
from  a  nozzle  is  well  known,12  as  ;.s  the  velocity  distribution.14  Con¬ 
sider  first  the  low  Reynolds  number  case. 

When  the  modified  Reynolds  number  is  less  than  20,  where 


and  is  the  nozzle  throat  length,  the  flow  in  a  nozzle  is  fully 

developed.  For  aspect  ratios  near  1,  the  velocity  profile  is  near 
parabolic,  especially  in  the  narrow  dimension.  When  there  is  no  effec¬ 
tive  potential  core,  the  parabolic  profile  very  closely  approximates  the 
two-dimensional  jet  profile;  hence,  one  may  expect  the  parabolic  profile 
between  the  planes  to  remain  parabolic,  especially  because  there  is 
little  or  no  pressure  gradient. 

The  loss  coefficient  is  defined  as  the  ratio  of  the  momentum 
reaching  the  splitter,  JXSp,  to  the  momentun  at  the  nozzle  exit,  Jg,  so 


12F.  Manion  and  T.  Drzewiecki ,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I_,  Harry  Diamond  Laboratories  (October  1974), 

14T.  M.  Drzewiecki,  Interpretation  of  Surface  Static  Pressure  Distri¬ 
butions  in  Fluid  Amplifier  Applications,  Harry  Diamond  Laboratories 
HDL-TR-1627  (July  1973). 


CXSP  =  JXSP^Js'  w^ere  t'ie  momentum  flux  of  a  parabolic  profile  in 
general  is  J  =  (6/15)bghpu£ .  The  ratio  between  two  parabolic  profiles 
is  then  merely  the  ratio  of  the  squares  of  center-line  velocity.  From 
equation  (8),  the  ratio  of  the  velocities  is 


u  / 
cf» 


Uc(x  =  0) 


Hence,  the  loss  coefficient  becomes 


JXSP 
CXSP  =  _J  — 


(47) 


(48) 


When  there  exists  a  potential  core,  the  net  loss  coefficient  is 
a  product  of  the  losses  in  the  core  region  and  the  decaying  region. 
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The  loss 
in  equation  (48), 


coefficient  in  the  decaying  jet  portion  is  the  same  as 
except  that  it  acts  over  a  shorter  distance. 
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A  similar  analysis  is  performed  for  the  momentum  loss  in  the 
potential  core  region,  except  that  now  the  velocity  profiles  are  not  of 
the  same  form.  However,  the  momentum  in  the  vertical  velocity  profile 
is  known  when  the  boundary  layer  displacement  and  momentum  thicknesses 
are  known.  The  momentum  flux  in  the  profile  per  unit  width  is 
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The  ratio  of  these  momenta  between  any  two  points  is  therefore 


J 


5 *  8. 

b  *  b 


s/  X-X2 


2/6 l_  _  11  .  JL\ 

0  \bs  bs  bs/  X=X 


For  a  zero  pressure  gradient  and  flow  over  a  flat  plate , 
momentum  integral  analysis15  tells  us  that 


t-ter- 


(52a) 


(52b) 


2.61  . 


When  these  values  are  substituted  into  equation  (51)  and  evaluated  at 
the  core  length  and  the  nozzle  exit,  the  result  is 


Xn+Xcore 
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where  Xn  is  the  effective  length  of  the  nozzle  if  it  were  a  flat  plate 
as  given  in  equation  (16). 

In  general,  however,  one  may  write  equation  (51)  as 


£(■ 


6  ‘  5/ 
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15H.  Schlichting ,  Boundary  Layer  Theory,  McGraw-Hill  Book  Co.,  New 
York  (1960). 
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The  quantity  [1  -  (6*/6)  -  (9/6)]  increases  in  a  favorable 
pressure  gradient  and  decreases  in  an  unfavorable  one.  Ft>r  small 
pressure  gradients,  the  boundary  layer  thickness  is  roughly  constant. 
The  pressure  gradient  at  the  nozzle  exit  is  always  favorable,  while 
downstream  one  would  expect  a  slightly  adverse  gradient.  Although  there 
are  no  comprehensive  data  on  the  gradients  t°  verify  equation  (55),  the 
form  of  this  loss  equation  is  important.  The  values  of  [1  -  (6*/6)  - 
(0/5)]  can  be  adjusted  to  give  correlation  of  the  final  results  with 
experimental  data  on  pressure  recovery.  Comparison  with  a  large  number 
of  different  cases  has  indicated  that,  for  aspect  ratios  greater  than  1, 
there  is  one  value  and,  for  aspect  ratios  less  than  1,  there  is  another 
as  follows: 


(56) 


where 


A  =  4. 0,  B  =  5.6,  o  _>  1 , 
A  =  B  =  3.0,  o  <  1. 


The  value  of  X  can  be  the  core  length  or  the  splitter  distance 
so  that  the  loss  coefficient  for  the  case  where  the  jet  has  a  potential 
core  for  its  entire  length  is 


c  = 

XSP 


(57) 


and  the  loss  coefficient  when  a  core  is  present  is 


(58) 
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The  validity  of  these  loss  coefficients  is  demonstrated  espe¬ 
cially  well  when  the  pressure  recovery  of  a  laminar  jet  device  is 
measured,  since  pressure  recovery  is  a  direct  measure  of  that  coeffi¬ 
cient  as  seen  from  equation  (24),  in  which 
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where 


C,  =  c.C_C 
loss  0  XSP  Bsp 


(60) 


Figures  8  to  10  illustrate  the  kind  of  agreement  obtained  with  the 
theory  and  experimental  data  for  various  splitter  distances,  aspect 
ratios,  and  output  widths.  These  equations  and  those  yet  to  be  devel¬ 
oped  have  been  incorporated  into  the  numerical  algorithm  proposed  by 
Drzewiecki . The  theoretical  curves  presented  in  figures  8  and  9  were 
generated  by  that  program. 


SUPPLY  PRESSURE  (mm  Hg) 

Figure  8.  Pressure  recovery  of  long,  xth  =  40,  nozzle 
laminar  jet  angular  rate  sensor:  normalized 
splitter  length  =  20,  supply  width  =  0.889  mm, 
normalized  output  width  =  1.7,  and  aspect 
ratio  =  1.1. 


16T.  M.  Drzewiecki,  Fluerics  38:  A  Computer  Aided  Design  Analysis  for 
the  Static  and  Dynamic  Port  Characteristics  of  Laminar  Proportional 
Amplifiers ,  Harry  Diamond  Laboratories  HDL-TR-1758  (June  1976). 
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Figure  9. 


I  100  200  300  400  500 

REYNOLDS  NUMBER 

Pressure  recovery  of  short  nozzle  laminar  jet  angular 
rate  sensor:  normalized  splitter  length  =  11.11,  supply 
width  =  1.14  mm,  normalized  output  width  =  1.67,  and 
aspect  ratio  =  2.77. 
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Figure  10.  Pressure  recovery  for  laminar  proportional  amplifiers  of 
different  aspect  ratios  (o's):  normalized  splitter 
length  =  8  and  supply  width  =  0.5  mm. 
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An  updated  version  of  the  program  is  given  in  appendix  A.  This 
program  is  a  general  laminar  jet  device  algorithm  and  generates  far  more 
information  than  the  rate  sensitivity.  The  input  and  output  character¬ 
istics  also  are  available.  This  program  now  gives  the  rate  sensitivity 
of  amplifiers  and,  conversely,  the  gain  of  rate  sensors. 

Figures  11  and  12  allow  comparison  between  theory  and  experi¬ 
ment  for  the  overall  sensitivity  of  an  LJARS.  The  agreement  is  con¬ 
sidered  acceptable. 


0  30  60  90  120  150  180  210 

MODIFIED  REYNOLDS  NUMBER 


Figure  11.  Sensitivity  of  laminar  jet  angular  rate  sensor  for 

various  configurations:  Xth  =  nozzle  throat  length, 
Xsp  =  splitter  length,  bs  =  supply  width,  and 
o  =  aspect  ratio. 


REYNOLDS  NUMBER 

Figure  12.  Sensitivity  of  laminar  jet  angular  rate 
sensor  for  hydraulic  operation. 

For  engineering  purposes,  the  simple  relation  for  the  losses 
developed  by  Manion  and  Drzewiecki12  for  >  50  may  be  used  without 
loss  of  accuracy.  The  expression  for  sensitivity  is  therefore 

12F.  Manion  and  T.  Drzewiecki,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I_,  Harry  Diamond  Laboratories  (October  1974). 
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A  listing  of  this  equation  is  contained  in  appendix  A. 


(61) 


2.4  Dynamic  Response 

The  response  of  the  jet  to  am  applied  sinusoidal  angular  rate 
is  complex.  Ideally,  with  no  field  effects,  the  jet  response  would  be 
governed  only  by  the  transport  delay.  In  other  words,  the  transfer 
function  would  be  a  pure  delay.  Such  is  not  the  case  in  reality.  It 
takes  the  field  a  finite  amount  of  time  to  respond  to  the  moving  of  the 
bounding  surfaces  so  that  the  device  reacts  like  an  ideal  sensor  for  the 
first  instant  of  time  when  subjected  to  a  step  input  of  rate,  but  then 
the  transient  of  the  field  catches  up.  Figure  13  shows  the  output 
differential  pressure  response  of  an  I^TARS  to  a  step  in  rate .  There  are 
regions  that  do  not  behave  as  one  may  have  expected.  For  example,  an 
initial  negative  spike  is  followed  by  an  increasing  signal ,  then  a 
second  but  slower  negative  region,  a  second-order  rise  to  a  slight 
overshoot,  and  a  settling  down  to  a  constant  value.  The  same  happens 
when  the  rate  is  removed. 

The  initial  negative  spike  is  attributable  directly  to  the 
sudden  acceleration  of  the  splitter,  which  protrudes  into  the  vent  area 
and  acts  like  a  paddle.  Figure  14  illustrates  that,  when  the  splitter 
moves  upward,  a  stagnation  (positive)  pressure  develops  on  the  top  side 
and  a  suction  (negative)  pressure  develops  below.  This  is  a  negative 
differential  pressure  transmitted  at  acoustic  speed  to  the  output  ports 
as  opposed  to  a  final  positive  pressure  differential  at  a  steady  rate  in 
which  the  jet  is  lagging  the  stru  ture  and  tends  to  develop  a  larger 
pressure  on  the  lower  output  and  a  smaller  pressure  at  the  upper  out¬ 
put.  The  paddle  pressure  signal  is  transmitted  at  the  speed  of  sound  to 
the  output  transducer . 

After  this  spike,  the  jet  should  lag  the  splitter  and,  hence, 
the  pressure  in  figure  13  rises.  To  explain  the  occurrence  of  the 
second  negative  part  requires  some  analysis.  However,  we  present  only 
some  phenomenological  observations  obtained  from  flow  visualization 
research.  Figures  15  and  16  show  the  instantaneous  jet  position  as 
marked  by  a  dye  stream  for  two  intermediate  times.  In  figure  15,  the 

time  corresponds  to  a  short  time  during  the  initial  rise  of  the  output 

pressure  in  figure  12.  In  figure  16,  the  time  is  at  an  instant  just 

about  at  the  minimum  of  the  second  negative  portion.  Qualitatively,  we 

observe  that  the  input  disturbance  at  the  nozzle  exit  travels  at  half 
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the  average  jet  velocity.  This  speed  is  apparent  when  one  considers  a 
disturbance  as  an  eddy  traveling  downstream  like  a  wheel  on  the 
ground.  The  part  of  the  eddy  in  contact  with  the  jet  has  the  axial 
velocity  of  the  jet,  yet  the  part  in  the  ambient  field  has  zero  veloc¬ 
ity;  hence,  the  centroid  of  the  eddy  must  travel  half  the  jet  veloc¬ 
ity.  Figure  17  illustrates  this  phenomenon. 


RATE  TABLE  TACHOMETER 


UARS  RESPONSE 


INPUT  RATE  STEP:  0  TO  +10  DEG/s 
SUPPLY  FLOW:’  1.89  /  10-5  m3/i 
TIMESCALE:  O.U/cm 


RATE  TABLE  TACHOMETER 


- p - 

NO  SUPPLY  FLOW 
INPUT  STEP:  +10DEG/I 

Figure  13.  Oscilloscope  traces  of  laminar  jet  angular 
rate  sensor  (LJARS)  step  response  in 
hydraulic  oil. 


Figure  14.  Paddle  effect  due  to  sudden  acceleration. 
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The  signal  arrives  at  the  splitter  in  twice  the  particle 
transport  time.  The  disturbance  as  it  travels  to  the  splitter  grows  and 
causes  the  jet  center  line  to  deflect  whiplike,  causing  the  negative 
part  of  the  output  response.  The  jet  settles  out  to  its  steady-state 
deflection  in  yet  another  disturbance  transport  time.*  A  steady  state 
is  achieved  in  two  disturbance  transport  times  or  four  average- 
particle-velocity  transport  times,  as  is  evident  from  figure  12.  We 
have  not  touched  on  the  reason  for  the  whiplike  action  of  the  jet.  This 
is  a  topic  for  further  research;  however,  we  provide  some  rationale  and 
direction  to  such  research  at  the  end  of  this  section. 

Ftar  design  engineering,  a  delay  of  four  times  the  transport 
time  predicts  a  phase  shift  of  -90  deg  at  a  frequency  far  below  that  of 
the  dynamic  effects  just  stated.  Since  phase  shift,  <f>,  is  defined  as 
the  ratio  of  delay  time,  Td,  to  period  of  oscillation,  then  in  degrees 

Td 

(j>  =  360  y  ,  (62) 

where  Td  is  four  times  the  average  transport  time,  and  X  is  the  period 
of  the  frequency  =  1/f,  so  that 


$  =  360f  T 

a 

and,  at  =  90  deg, 


(63) 


(64) 


*The  first  edge-tone  frequency  occurs  at  a  frequency  roughly 
corresponding  to  twice  the  transport  time.  If,  indeed,  a  disturbance 
travels  to  the  splitter  at  that  speed  and  is  acoustically  (virtually 
instantaneously )  fed  back,  then  it  is  reasonable  that  the  first 
edge-tone  frequency  should  be  at  0.466Vs/XSp,  where  Vs  is  the  supply 
velocity.  Higher  harmonics  do  not  occur  as  multiples  because  at  higher 
frequencies  the  acoustic  delay  becomes  appreciable  and  changes  the  mode. 
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Remembering  that 


N 


R 


V 


1/2 


we  have 


f 


<j>=-90# 


c  Nv 
d  R 


16X 


b2 
sp  s 


(65) 


From  this  argument,  we  may  conclude  that  the  frequency  response  of  an 
LJARS  is  flat  past  the  -90  deg  phase  shift  point.  Experimental  evidence 
bears  this  out  quite  adequately  as  shown  in  figures  18  to  20.  Fig¬ 
ures  18  and  19  show  the  Bode  response  of  an  UARS  used  in  a  tank  gun 
stabilization  system  operating  in  hydraulic  oil.  It  can  be  shown  that 
the  phase  shift  in  figure  18  as  determined  from  equation  (65)  is  in 
excellent  agreement,  and  consequently  the  point  at  which  -90  deg  occurs 
is  well  defined  by  the  above  theoretical  considerations.  Examination  of 
equation  (65)  indicates  that  devices  of  the  same  size  operating  at  the 
same  Reynolds  number  have  the  same  actual  bandwidth,  regardless  of  fluid 
medium,  provided  that  the  kinematic  viscosity  remains  the  same.  At  room 
temperature,  the  kinematic  viscosities  of  MIL-H-5606  hydraulic  fluid  and 
air  are  the  same.  So  a  hydraulic  rate  sensor  and  a  pneumatic  rate 
sensor  of  the  same  size  have  the  same  response  and  bandwidth.  This  is  a 
different  effect  than  was  found  for  sensitivity  where  sensitivity  is 
shown  to  be  independent  of  size. 


Figure  18.  Dynamic  response  of  hydraulic  laminar 

jet  angular  rate  sensor  (from  AiResearch) . 
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SUPPLY  FLOW  =  0.18 gpm 


Figure  19.  Dynamic  response  of  hydraulic  laminar  jet 
angular  rate  sensor. 
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20.  Dynamic  response  of  overscale  4:1  hydraulic  laminar 
jet  angular  rate  sensor. 


For  the  whiplash  jet  response,  if  it  is  assumed  that  the 
acceleration  effects  settle  out  before  the  rate  sensor  can  respond,  then 
the  fluid  in  the  sensor  field  takes  on  the  solid  body  rotation  of  the 
sensor  body  with  appropriate  secondary  motion  as  required  by  the  govern¬ 
ing  equations.  For  no  acceleration  effect,  the  jet  stream  enters  the 
field  with  a  velocity  relative  to  the  field.  This  velocity  is  independ¬ 
ent  of  the  center  of  rotation  of  the  sensor  body.  Experiments  performed 
on  an  over scale  hydraulic  sensor  gave  the  same  response  trace  as  shown 
in  figures  13  and  15,  whether  the  sensor  was  pivoted  at  the  nozzle  exit 
plane  at  the  splitter  or  half  way  between  the  nozzle  exit  and  the 
splitter.  These  results  support  the  assumption  that  acceleration 
effects  settle  out  before  the  sensor  responds  to  rate  input. 

For  this  assumption,  the  whiplash  response  shown  in  figures  15 
and  16  can  be  qualitatively  examined.  As  noted  in  figure  3,  the  jet 
particles  enter  the  field  with  a  relative  velocity  and  tend  to  move  out 
of  the  former  jet  path  shown  in  figure  21  by  the  dashed  lines.  The 
nozzle  and  the  splitter  also  have  moved  in  the  time  increment  since  the 
rate  was  applied. 
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Figure  21.  Transient  jet  particle  path  under  step  rate  response; 

motion  is  frozen  at  time  increment  after  step  in  rate 
was  applied. 


The  jet  stream  is  still  continuous  and  includes  particles  that 
had  left  the  nozzle  before  the  rate  was  initiated.  The  resulting  jet  is 
as  sketched  in  figure  21 ,  and  it  gives  rise  to  the  formation  of  two 
large  eddies.  The  leading  eddy  is  formed  as  a  result  of  the  displace¬ 
ment  of  the  jet  in  the  field  and  is  essentially  a  startup  vortex.  Ttoe 
law  of  conservation  of  angular  momentum  requires  that  the  diameter  of 
the  eddy  increase  with  decreasing  edge  velocity.  Since  the  jet  edge 
velocity  decays  downstream,  it  follows  that  the  eddy  increases  in 
diameter.  The  starting  vortex  is  formed  on  the  opposite  side  of  the  jet 
from  the  signal  eddy  and,  hence,  it  draws  the  jet  particles  across  the 
instantaneous  jet  center  line.  This  then  accounts  for  the  negative 
signal  when  that  eddy  reaches  the  splitter. 

Considerable  analysis  is  still  required  of  both  the  amplitude 
dynamics  and  the  phase  dynamics  of  the  LJARS.  Correlation  of  step 
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response  analysis  to  sinusoidal  response  predictions  also  is  neces¬ 
sary.  What  is  presented  here  is  merely  an  engineering  guide  to  design. 


3.  SENSOR  DESIGN  CONSIDERATIONS 

This  section  deals  with  the  analyses,  the  design  guides,  and  the 
rules  of  thumb  necessary  to  implement  an  LJARS  in  a  system. 

3.1  Operating  Range 

The  operating  range  of  the  LJARS  is  determined  by  two  limits. 
Primarily,  the  useful  upper  limit  occurs  at  the  onset  of  turbulence  in 
the  field.  This  has  been  established^2  for  devices  having  fairly  short 
nozzle-to-splitter  distances  (8  £  xgp  _£  10).  The  basic  phenomenon  of 
transistion  to  turbulence  is  complex  and  is  buried  in  empiricism. 
However,  for  laminar  jet  devices,  the  phenomenon  should  follow  the 
relation  for  the  onset  of  instability  where  there  exists  a  critical 
Reynolds  number  based  on  displacement  thickness  for  each  specified 
pressure  gradient.  Schlichting1 5  discusses  this  at  length,  as  does 
Hinze.*®  In  principle,  however,  for  low  or  almost  zero  pressure 
gradients,  the  result  is 
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This  value  may  be  as  low  as  420.  From  equation  (52b), 


We  may  thus  establish  a  critical  nozzle  Reynolds  number.  X  is 
the  effective  flat  plate  distance  to  the  splitter,  X  =  Xgp  +  xeff 
(effective  length).  xeff  is  obtained  from  equations  (15)  and  (16). 
This  results  in 


XN„  ...  , 

R  critical 


1.4  x  105  . 


(67) 


12F.  Manion  and  T.  Drzewiecki,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I_,  Harry  Diamond  laboratories  (October  1974), 

15H.  Schlichting ,  Boundary  Layer  Theory,  McGraw-Hill  Book  Co.,  New 
York  (1960). 

l8J.  0.  Hinze,  Turbulence,  McGraw-Hill  Book  Co.,  New  York  (1959). 
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For  Xflp  ■  8  and  Xejf  *  2,  X  -  10?  this  leads  to  a  critical  Reynolds 
number  of  14,000,  which  is  about  10  times  too  high.  The  adverse 
pressure  gradient  at  the  splitter  should  cause  earlier  transition.  As 
we  note  from  Schlichting, 18  the  critical  Reynolds  number  may  easily  be 
an  order  of  magnitude  lower  in  an  adverse  pressure  gradient.  Experi¬ 
ments  have  shown  consistently  that,  for  laminar  devices,  the  relation¬ 
ship 


°^R  critical 


1400 


(68) 


holds  true  with  the  condition  that  the  flow  entering  the  devices  be  free 
from  noise  or  turbulence.  (This  requires  a  fairly  smooth  manifold, 
sect.  3.7.)  For  a  -  1,  this  would  indicate  that 


^R  critical 


1.4  x  104 


or 


ON 


R  critical 


=  1.4  x 


104  S 


(69) 


a  can  be  replaced  with  an  observed  value  for  6*/bg  -  0.06o1j/2  for 
0.2  <  o  <  4  with  equations  (52b)  and  (69),  which  yields 

crltlc.1  '  *  s(’°6)]  1/2  -  25»V  •  (’0) 

This  is  shown19  plotted  in  figure  22.  Now  Hinze18  states  that,  for 
Oouette  flow  between  parallel  plates,  the  criterion  for  transition  is 


uh 

v 


1100  . 


(71) 


Noting  that  u  =  NR  =  umaxb  g/v  (where  U  is  the  mean  or  reference 

velocity),  and  o  =  h/b8. 


15tf.  Schlichting ,  Boundary  Layer  Theory,  McGraw-Hill  Book  Co.,  New 
York  (1960). 

18J.  0.  Hinze,  Turbulence ,  McGraw-Hill  Book  Co.,  New  York  (1959). 

19G.  Roe,  Fluidic  Laminar  Angular  Rate  Sensor  Research  Program, 
McDonnell  Douglas,  Titusville ,  FL,  HDL-CR-7 5-020-1  (November  1975). 
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Figure  22.  Critical  Reynolds  number  as  function  of 
splitter  distance. 


At  the  point  of  transition,  the  discharge  coefficient  of  short 
nozzles  is  on  the  order  of  0.75  to  0.8,  which  leads  to 

1375  <  oN  <  1470  ,  (72) 

R 


which  also  is  shown  in  figure  22.  The  flow  must  be  stable  if  the 
Reynolds  number  is  below  the  Couette  flow  critical  number  since  all  the 
flow  is  at  or  below  that  NR.  Hence,  the  flows  that  fall  above  both 
curves  are  transitional  (or  turbulent).  For  the  range  of  interest,  X 
>_  8,  we  may  use  the  simple  expression  in  equation  (68).  The  rule  of 
thumb  that  is  used  in  the  laboratory  is 


R  transition 


1400  . 
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There  is  no  effect  of  splitter  distance  until  the  splitter  distance  is 
less  them  8;  then  the  flow  is  more  stable,  as  one  might  expect.  The 
scatter  of  the  data  is  due  to  aspect  ratio  variations.  Recent  experi¬ 
mental  evidence  correlates  aspect  ratio  and  the  theoretical  arguments  of 
Manion  and  Drzewiecki .*2  This  leads  to  a  more  universal  parameter  for 
transition  wherein 


Just  before  transition  when  the  offset  is  bad,  an  LJARS  tends 
to  exhibit  digital  or  erratic,  albeit  laminar,  behavior.  This  clearly 
limits  the  useful  range  of  such  a  device.  Figure  23  shows  a  plot  of 
differential  output  as  a  function  of  supply  pressure  for  a  device  with 
poor  offset  characteristics.  It  is  clear  that  the  onset  of  turbulence 
occurs  at  about  14  in.  of  H20;  however,  it  is  undesirable  to  operate 
this  device  in  the  region  where  the  change  of  output  with  supply  pres¬ 
sure  is  so  drastic  and  erratic.  It  is  far  more  advantageous  to  operate 
on  the  flat  portion  between  8  and  10  in.  of  H20.  This  obviously  means 
that  the  sensitivity  will  be  lower  at  lower  NR.  Experimental  data  indi¬ 
cate  that  a  30-percent  backoff  from  critical  NR  normally  puts  one  in  the 
range  of  acceptable  performance.  As  shown  in  figure  23,  a  reduction  of 
supply  pressure  by  a  factor  of  two  from  the  transition  point  of  14  in. 
of  H20  (a  30-percent  reduction  in  NR)  puts  the  device  at  7  in.  of  H20, 
which  appears  to  be  a  much  better  place. 


For  safe  operation,  therefore,  we  suggest  a  maxinmn  operating 
Reynolds  number  defined  by 


oNr  =  1000  .  ,  (73) 

The  lower  limit  must  necessarily  be  defined  by  the  user  at  the 
point  below  which  the  sensitivity  is  too  low  for  use.  As  a  rule,  a 
50-percent  reduction  of  sensitivity  may  be  acceptable.  Figure  11  shows 
typical  sensitivity  as  a  function  of  Reynolds  number  for  various 
devices.  If  the  requirement  for  only  50-percent  change  in  sensitivity 
is  imposed,  the  device  depicted  by  the  left-most  data  may  operate 
between  NR  =  8  and  NR  *  7  or  a  supply  pressure  change  of  only  a  factor 
of  1.3.  This  operation  is  caused  primarily  by  the  very  long  nozzle 
length,  whose  flow  characteristics  are  extremely  sensitive  to  change  in 
Nr.  On  the  other  hand,  a  properly  designed  short  nozzle  configuration. 


12f.  Manion  and  T.  Drzewiecki ,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers ,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  I_,  Harry  Diamond  Laboratories  (October  1974). 
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while  giving  a  similar  rate  sensitivity,  demonstrates  considerably  less 
NR  sensitivity.  Such  a  device  may  operate  from  as  high  as  =  120  to 
as  low  as  NR  =  75,  which  is  a  factor  of  2.56  in  pressure  or  a  good 
factor  of  5  better  in  percentage  than  the  latter  case  (156  versus  30 
percent) .  This  means  that  supply  pressure  regulation  may  not  have  to  be 
as  stringent  a  requirement. 


Figure  23.  Null  offset  for  early  laminar  jet  angular  rate  sensor 
(from  McDonnell  Douglas  (19)). 


For  engineering,  we  may  allow  NR  to  vary  by  a  factor  of  two 
(pressure,  by  a  factor  of  four).  (This  incidentally  corresponds  to  the 
Reynolds  number  change  imposed  when  operating  a  pneumatic  system  between 
-55  and  +80  C  or  -65  and  +180  F.)  We  may  thus  arbitrarily  define  the 
LJARS  operating  range  as 


500  <  aN  <  1000  , 

K 


(74) 


and  this  is  the  range  referred  to 


3.2  Environmental  Sensitivity  and  Compensation 


Environmental  sensitivity  concerns  itself  with  changes  in 
operating  characteristics  due  to  changes  in  the  environment.  The 
environment  consists  of  many  different  parameters  such  as  temperature, 
pressure,  contamination,  noise,  and  other  perturbations;  this  section 
considers  only  those  effects  that  cause  changes  in  the  operating 
Reynolds  number — temperature  and  pressure.  The  other  effects  are 

considered  in  later  sections. 

We  briefly  mention  temperature  effects  in  section  3. 1  during 
the  discussion  of  operating  range.  Temperature  and  pressure  (altitude) 
affect  the  Reynolds  number,  Np  =  pbgUg/u,  through  their  influence  on 
both  density  and  viscosity.  When  the  operating  pressure  or  flow  is  held 
constant  in  an  application,  the  Reynolds  number  changes  due  to  the 
environmental  changes.  Viscosity  changes  with  temperature,  and  density 
changes  with  temperature  and  pressure  or  altitude.  The  predominant 
effects  are  two-fold;  (1)  the  sensitivity  changes,  and  (2)  the  oper¬ 
ating  point  moves  with  respect  to  the  point  of  transition  to  turbu¬ 
lence.  In  any  application,  the  limits  of  the  Reynolds  number  must  be 
determined,  and  the  operating  point  must  be  set  so  that  it  falls  within 
an  acceptable  range. 

Changes  in  sensitivity  are  crucial  in  maintaining  the  stability 
of  a  feedback  control  system .  Too  much  gain  often  is  worse  than  not 
enough  because  the  system  can  become  unstable.  Low  gain  prevents  the 
system  from  functioning.  Transition  to  turbulence  increases  the  system 
noise,  degrades  the  gain,  and  increases  the  threshold  levels,  perhaps  to 
an  undesirable  point.  Typically,  the  signal-to-noise  ratio  for  turbu¬ 
lent  devices  may  be  as  low  as  50:1  as  opposed  to  106:1  for  laminar 
devices. 


To  compensate  an  LJARS  for  temperature,  one  must  either  regu¬ 
late  the  Reynolds  number  or  change  the  gain  outside  the  parameters  of 
the  LJARS  operation. 

Consider  regulation  of  the  Reynolds  number  when  viscosity  is 
changing.  When  constant  mass  flow,  m,  is  supplied  to  an  orifice,  m0<- 
and  a  capillary,  m  ,  in  parallel  as  shown  in  figure  24,  the  distri¬ 
bution  of  flow  depends  on  the  relative  fluid  viscosity  and  density  since 
the  resistance  of  the  orifice  depends  on  density  only,  whereas  that  of 
the  capillary  depends  on  viscosity  only.  As  the  viscosity  increases 
(the  temperature  decreases  in  a  liquid  or  increases  in  a  gas),  less  flow 
passes  through  the  capillary  and,  hence,  more  goes  through  the  ori¬ 
fice.  The  equation  for  the  conservation  of  mass  is  written  for  the 
resistor  pair : 
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so  that 


or 
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where 


R’  =  AP/m  . 


P,  m 


PRESSURE,  MASS 
FLOW  RATE 


Figure  24.  Parallel  orifice  (o)- 
capillary  (ca)  pair. 


Now  the  orifice  resistance  depends  on  flow,  R  =  K  m  ,  orifice  coeffi¬ 
cient  KQ  =  (1/2)c^A2p,  and  A  is  the  area,  so  that  °  ° 


m_ 


PQ  = 


Rc 


(76 


As  the  viscosity  (or  R  )  increases,  so  does  m  . 

c  o 

The  ratio  of  the  flow  at  one  viscosity  to  that  at  another  vis¬ 
cosity  is 
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(77) 


In  the  limit  as  the  capillary  resistance  becomes  very  small  (one  by¬ 
passes  a  large  amount  of  flow  around  the  orifice) ,  R  <<  K  m, 

c  o 


*02 


1/2 


(78) 


Capillary  resistance  is  directly  related  to  viscosity;  hence. 


Now  the  flow  Reynolds  number  is  defined  as  N  =  m/hu;  hence,  the  ratio 
of  two  orifice  Reynolds  numbers  is 


N 


R2 


o2 


(80) 


Substituting  equation  (79)  into  equation  (80),  therefore,  yields 


*R1 


"R2 


1/2 


(81) 


This  result  clearly  indicates  that,  as  viscosity  increases,  the  Reynolds 
number  still  decreases;  however,  instead  of  decreasing  linearly  with  the 
reciprocal  of  viscosity,  it  now  decreases  as  the  square  root  of  that 
quantity.  For  a  change  in  viscosity  of  a  factor  of  four,  the  Reynolds 
number  changes  only  a  factor  of  two  under  this  condition. 

One  cannot  pay  the  penalty  of  infinite  flow  bypass,  so  that 
equation  (82)  must  be  used. 
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(82) 


We  now  recall  equation  (26)  for  the  rate  sensitivity  of  the  LJARS  that 
has  an  orifice-like  supply  nozzle. 


2  NR  Prec  _ 

S  =  573  57  —  ZOy  ' 

d  s 


and  we  observe  the  ratio  of  two  sensitivities  at  two  viscosities. 


(83) 


If  we  linearly  approximate  the  pressure  recovery  in  terms  of  Reynolds 
number  (fig.  9)  so  that  Prec/Ps  =  kNR,  where  k  is  a  constant  and  assume 
that  Z01  =  ZOj ,  then 


"k  12 

u2 

R2 


(84) 


no  compensation 


In  the  best  compensated  case,  equation  (81),  we  substitute  for  N  /N 

(noting  that  N  =  c  N  )  R 

R  u  R 


Ibis  says  that,  as  the  Reynolds  number  increases  (c^  increases  with  NR)  , 
sensitivity  decreases.  Since  we  have  seen  that  a  temperature-related 
Reynolds  number  increase  cannot  be  halted  entirely,  then  sensitivity  can 
decrease.  Without  compensation,  the  sensitivity  increases  with  increas¬ 
ing  Reynolds  number  (increasing  c^)  .  When  one  considers,  however,  that 
the  supply  nozzle  of  a  sensor  is  not  a  pure  orifice,  that  a  shunt  capil¬ 
lary  is  not  a  pure  linear  resistor,  and  that  one  cannot  afford  to  shunt 
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an  infinite  amount  of  flow,  then  the  aatount  of  compensation  is  consider¬ 
ably  reduced.  This  reduced  compensation  can  be  determined  by  solving 
the  complete  equation  for  the  flow,  equation  (76).  One  can  rewrite 
equation  (76)  for  the  Reynolds  number,  Nr  *  m0/c,jhy,  and  define  ratios 
of  viscosity,  resistance,  and  discharge  coefficient  relative  to  some 
nominal  conditions  identified  by  subscript  "o."  Also,  the  total  flow  is 
held  constant  so  that 


=  N_  c.  u  h 
Ro  do  o  o 


where 


R  1  =  K  m  I 

o|o  o  o Jo 

We  further  define  the  ratio  of  nozzle  resistance  to  shunt 
resistance  at  the  nominal  point  R^/ R^,  I  Q .  After  considerable  manipu¬ 
lation,  equation  (76)  becomes 


(85) 


The  universal  equation  for  the  discharge  coefficient  is  given  as  equa¬ 
tion  (86)  and  shown  in  figure  25, 


R  =  [6XP  (6*6Cd)  '  ’]  ' 


(86) 


so  that  these  two  equations  may  be  successfully  iterated.  An  iterative 
computer  program  is  given  in  appendix  A.  The  result  is  perfectly  gen¬ 
eral  for  all  linear  shunts  (capillaries).  Figure  26  shows  the  results 
generated  from  equations  (85)  and  (86)  for  a  nozzle  with  a  high  dis¬ 
charge  coefficient,  meaning  that  the  linear  portion  is  small.  The 
efficacy  of  such  a  scheme  can  be  observed.  However,  this  compensation 
scheme  is  of  no  value  when  the  nominal  operating  modified  Reynolds 
number  is  low  for  a  very  long  nozzle,  such  as  40  nozzle  lengths.  Fig¬ 
ure  27  shows  the  nozzle  Reynolds  number  as  a  function  of  shunt  resist¬ 
ance,  illustrating  the  futility  of  shunting  in  such  a  case.  Even  on  an 
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expanded  scale,  little  change  Is  observed.  (Merely  using  a  constant 
flow  source  in  itself  compensates  a  little  for  viscosity  changes.  If 
pressure  is  held  constant,  when  viscosity  changes  by  a  factor  of  five, 
so  does  the  Reynolds  number;  however,  at  a  constant  flow,  the  same 
factor  of  five  yields  only  a  factor  of  four,  as  shown  in  figure  26.  It 
does  so  due  to  a  slight  increase  in  supply  pressure  with  increasing 
viscosity  for  constant  flow  due  to  an  increased  overall  resistance.. 


Figure  25.  Universal  discharge  coef¬ 
ficient  (c^)  function. 


Figure  28  shows  how  kinematic  viscosity  changes  for  air  and  two 
types  of  hydraulic  fluid.  Hydraulic  fluid  dramatically  changes  in  vis¬ 
cosity  over  a  typical  vehicular  temperature  range  (6:1  to  17:1),  whereas 
air  changes  less  than  a  factor  of  two  over  the  same  range.  When  the 
rule  of  thumb  is  applied  that  the  Reynolds  number  excursion  can  be  kept 
to  half  the  viscosity  excursion,  hydraulic  fluid  may  still  change  from 
3:1  to  8:1.  In  this  manner,  by  compensating  for  viscosity  changes,  we 
compensate  for  changes  in  temperature. 
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VISCOSITY  (/>0/p) 


Figure  26.  Change  of  Reynolds  number  (NR/NRo)  with 
viscosity  (p0/p)  for  different  shunt 
resistances . 


In  figure  29,  data  show  the  amount  of  compensation  possible 
with  a  shunt  resistor  roughly  equal  to  seven  times  the  room  temperature 
nozzle  resistance.  When  the  flow  and  the  temperature  are  high,  the 
sensitivity  can  be  reversed  as  predicted.  In  this  case,  at  about  47  C, 
the  shunt  resistance  is  roughly  equal  to  the  nozzle  resistance  so  that, 
above  that  temperature,  reversal  compensation  occurs. 


When  viscosity  decreases,  LJARS  sensitivity  may  be  made  to 
decrease.  Also,  among  the  typical  gain  characteristics  of  fluidic 
amplifiers,  the  gain  increases  with  decreasing  viscosity  (an  increasing 
modified  Reynolds  number)  as  shown  in  figure  30.  It  is  therefore  pos¬ 
sible  to  achieve  a  flat  or  compensated  net  sensor-gain-block  rate  sensi¬ 
tivity  over  a  wide  band  of  temperatures. 


VISCOSITY 


Figure  27.  Effect  of  shunt  temperature  compensation 
on  long  nozzle. 


In  figure  31,  a  composite  gain  or  sensitivity  graph  is  deter¬ 
mined  from  the  gain  data  of  figure  30  and  compensated  sensitivity  from 
equation  (84)  and  figure  26.  The  sensitivity  of  a  pneumatic  rate  sensor 
is  superimposed  on  (multiplied  by)  the  gain  of  a  three-stage  gain 
block.  As  is  evident,  the  range  for  a  ±10  percent  variation  of  an 
uncompensated  system  would  be  unacceptable;  however,  when  the  LJARS  is 
compensated  by  using  shunt  flow,  the  system  can  operate  over  a  ±60  C 
range  (the  military  environment  range)  with  less  than  a  ±10  percent 
sensitivity  variation.  This  type  of  design  guide  can  be  generated  from 
the  following  rules  of  thumb: 

a.  The  gain  decreases  are  within  10  percent  per  stage  if  the 
Reynolds  number  decreases  less  than  50  percent  from  maximum  for  standard 
laminar  proportional  amplifiers  (LPA's)  (configuration  in  app  A).  (The 
gain  increases  10  percent  when  the  Reynolds  number  increases  100  percent 
from  midrange.) 


b.  Tfte  net  staged  gain  change  is  the  product  of  the  number 
of  stages  with  the  per- stage  gain  change. 

c.  The  Reynolds  number  changes  inversely  with  viscosity. 

d.  The  military  temperature  range  corresponds  to  a  factor  of 
two  in  viscosity  for  pneumatic  systems. 


•0  0  SO  160  260 

TEMPERATURE  (P) 


-*0  0  SO  100  ISO 

TEMPERATURE  (C) 

Figure  28.  Viscosity  versus 
temperature  for 
various  fluids. 

It  is  of  considerable  benefit  to  have  sensors  and  amplifiers 
that  exhibit  good  orifice  characteristics,  not  only  from  the  standpoint 
of  compensation,  but  also  because  they  inherently  are  less  sensitive  to 
temperature  change  and  have  higher  gain  and  better  overall  performance. 

In  conclusion  of  this  section  on  temperature  compensation, 
there  are  other  methods  to  control  the  Reynolds  number.  Drzewiecki20 
reports  a  bimetallic  resistive  element  designed  to  keep  the  Reynolds 

20T.  M .  Drzewiecki,  Fluidic  Resistive  Element,  U.S .  Patent  3,980,103 
(14  September  1976). 


number  constant.  Also,  classical  in  gaseous  systems  is  the  use  of  a 
parallel  orifice-capillary  pair.  The  orifice  resistance  depends  only  on 
density,  which  decreases  with  temperature,  and  the  capillary  resistance 
depends  on  viscosity,  which  increases  with  temperature.  In  such  a 
manner,  there  exists  a  region  of  temperatures  in  which  the  net  parallel 
resistance  is  constant.  Such  passive  compensation  componentry  is 
necessary  when  bridges  and  other  ancillary  circuitry  are  used  with  the 
LJARS  and  the  LPA.  Compensation  of  LPA's  may  be  achieved  by  uting 
operational  amplifiers  where  the  system  gain  depends  on  the  ratios  of 
feedback  to  input  resistors  (provided  that  there  is  sufficient  forward 
gain).  When  these  resistors  both  are  of  the  same  type,  there  is  no 
compensation.  When  the  feedback  resistor  depends  on  viscosity  and  the 
input  depends  on  density,  the  gain  varies  with  kinematic  viscosity 
(increases  with  temperature  in  air).  This  gain  may  be  shaped  by 
appropriate  selection  of  combinations  of  resistors.  Mon21  has  demon¬ 
strated  a  constant  gain  (to  ±10  percent)  from  modified  Reynolds  numbers 
of  20  to  above  100  or  a  range  of  5  to  1.  The  penalty  to  be  paid  here  is 
the  large  number  of  amplifiers  required  to  get  the  desired  gain  and, 
hence,  increased  complexity,  size,  and  flow.  Available  on  the  market 
are  temperature  compensated  flow  controllers.  Mechanical  devices  such 
as  pressure  and  flow  regulators  may  have  their  characteristics  adjusted 
to  provide  some  Reynolds  number  regulation.  Mechanical,  less  reliable 
devices  may  have  to  be  resorted  to  in  the  extreme  cases  of  hydraulic 
fluid  operation  of  from  -40  to  +80  C,  for  example.  In  such  an  event, 
the  system  reliability  would  be  limited  by  such  a  device. 

Changes  in  density  may  be  compensated  for  by  enclosing  the 
system,  pressurizing  it,  and  allowing  it  to  bleed  to  ambient  conditions 
through  a  sonic  orifice  as  suggested  by  the  work  at  McDonnell 
Douglas. 19,22  In  this  manner,  any  pressure  disturbances  outside  are  not 
able  to  propagate  into  the  system.  This  also  successfully  isolates  the 
system  from  outside  noise  in  addition  to  altitude  changes.  Care  must  be 
exercised,  however,  to  insure  that  the  sonic  orifice  does  not  itself 
generate  noise  due  to  turbulent  flow  upstream  of  the  throat. 


19G.  Roe,  Fluidic  Laminar  Angular  Rate  Sensor  Research  Program, 
McDonnell  Douglas,  Titusville,  FL,  HDL-CR-7 5-020-1  ( November  1975). 

2 George  Mon,  Flueric  Laminar  Gain  Blocks  and  an  Operational  Amplifier 
Scaler,  Harry  Diamond  Laboratories  HDL-TR-1730  (December  1975). 

22W.  J.  Westerman  and  R.  E.  Wright,  Evaluation  Program — Fluidic 
Laminar  Rate  Sensor — Final  Report,  McDonnell  Douglas,  Titusville,  FL, 
NAVAIR  L0252  (8  March  1974). 
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Figure  29.  Temperature  compensation  of  laminar  jet  angular 
rate  sensor  (LJARS) :  normalized  splitter 
length  =  20,  supply  width  =  0.75  mm,  aspect 
ratio  =  2.5,  and  shunt  resistance  =  7  *  LJARS 
supply  resistance. 
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Figure  31.  Temperature  compensation  by 
capillary  shunt  resistance 
of  subsystem  of  laminar  jet 
angular  rate  sensor  (LJARS) 
and  laminar  proportional 
amplifier  (LPA) . 


3.3  Shock  and  Vibration  Sensitivity 

As  a  follow  on  to  environmental  sensitivity,  it  is  of  interest 
to  consider  the  effects  of  shock  and  vibration. 

From  a  theoretical  standpoint,  shock  and  vibration  constitute 
accelerations  in  axes  perpendicular  and  parallel  to  the  laminar  jet. 
They -act  equally  on  all  fluid  particles  in  the  region.  As  long  as  there 
exist  no  density  gradients,  the  body  forces  produced  act  equally  on  all 
particles;  thus,  no  net  forces  act  on  the  jet.  In  addition,  those 
accelerations  along  the  jet  axis  and  those  perpendicular  to  the  bounding 
planes  produce  body  forces  that  are  in  the  common  mode  to  the  outputs. 
They  would  not  contribute  any  differential  signal  if  the  common  mode 
rejection  ratio  (CMRR)  were  infinite.  They  may  add  some  signal  in  the 
real  case;  however,  since  these  forces  do  not  move  the  jet,  the  CMRR  for 
such  a  phenomenon  is  extremely  high.  In  the  cross  axis,  there  exists 
the  possibility  of  jet  movement;  however,  since  the  structure,  the  field, 
and  the  jet  move  together  again,  no  differential  signals  should 
appear.  Rotational  accelerations  that  are  within  the  bandwidth  of  the 
device  will  be  detected.  In  section  2.4,  the  bandwidth  is  discussed. 
From  the  typical  frequency  response  of  figure  20  and  the  analysis,  we 
can  estimate  the  bandwidth  to  be  at  least  twice  the  -90  deg  of  phase 
shift  frequency,  probably  more  like  10  times.  Since  the  LJARS  measures 
rotational  rate,  however,  it  is  reasonable  that  rotational  shocks  and 
vibrations  should  pass  through.  Alternatively,  one  cam  assume  that 
noise  or  signals  come  from  angular  disturbances  only. 

As  further  demonstration  that  LJARS  characteristics  remain 
unchanged  in  vibrational  environments,  experiments  were  conducted  at 
AiResearch23  on  pneumatic  devices  in  the  frequency  range  of  5  to  1000  Hz 
at  an  amplitude  of  10  g,  in  all  three  orthogonal  axes.  These  experi¬ 
ments  showed  no  increase  in  noise  whatever  (as  monitored  on  oscillo¬ 
scopes)  and  no  effect  on  gain.  Analysis  and  experiments  on  an 
early  design  LJARS  at  General  Electric21*  coupled  with  laminar  amplifiers 
indicate  no  vibration  sensitivity  at  5  g  from  5  to  2000  Hz.  Shock 
loading  of  12,500  g  ( nonoperational)  did  not  alter  characteristics  of  or 
damage  the  LJARS.  The  most  susceptible  regions  are  internal  cavities 
that  are  excited  by  outside  perturbations.  In  general,  in  both  pneu¬ 
matic  and  hydraulic  operations,  these  resonances  are  at  very  high 
frequencies,  greater  than  5  kHz  for  devices  with  nozzle  dimensions  on 
the  order  of  1  ran. 

2 3Af.  Cycon,  T.  Sutton,  and  H.  Gamble,  Environmental  Test  Program  for 
the  Evaluation  of  Fluidic  Elements,  AiResearch  Manufacturing  Co.  of 
Arizona,  Phoenix,  AZ,  HDL-CR-75-187-1  (12  December  1975). 

2hv.  Bravo,  Feasibility  Investigation  of  a  Laminar  Rate  Sensor  (LARS), 
Phase  III  Report,  General  Electric  Co,  Schenectady,  NY,  71SD2118 
(December  1971). 
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3.4  Null  Offset  and  Null  Offset  Oompensation 


Perhaps  the  most  trying  aspect  of  fluidics  is  the  persistent 
presence  of  null  offset  as  shown  in  figure  23.  In  an  LJARS,  this  is  a 
false  rate  signal.  The  primary  causes  of  null  offset  are  geometric 
asymmetries.  If  the  jet  is  displaced  relative  to  the  outputs,  there 
exists  a  differential  output  pressure,  lhis  displacement  may  be  due  to 
a  misalignment  of  the  splitter  with  the  nozzle,  or  it  may  be  due  to  an 
angular  error  of  alignment  of  the  nozzle  axis  with  that  of  the 
splitter.  From  equation  (59),  it  is  apparent  also  that  the  recovered 
pressure  depends  on  the  output  width.  If  the  two  outputs  differ  in 
area,  then  a  pressure  difference  exists  with  no  jet  deflection.  Unequal 
control  resistance  or  edges  lead  to  a  differential  pressure  across  the 
jet,  which  causes  it  to  deflect.  This  deflection  depends  on  gain,  so 
that  at  low  Reynolds  numbers  there  is  little  deflection,  and  at  high 
Reynolds  numbers  there  is  large  deflection. 

Deflection  is  caused  also  by  asymmetric  separation  off  uneven 
nozzle  exit  corners  as  described  by  Manion  and  Drzewiecki. 12  Composite 
diagrams  of  these  effects  are  shown  in  figures  32  and  33.  One  may 
estimate  the  amount  of  offset  by  knowing  such  things  as  the  percentage 
of  the  difference  in  output  areas.  From  equation  (59),  it  follows  that 
there  is  a  one-to-one  relationship.  The  percentge  of  the  difference  in 
areas  gives  the  percentage  of  the  difference  in  recovered  pressures. 

It  is  not  so  simple  with  deflection,  however.  Consider  that 
for  an  LJARS  with  Xgp  =  20,  full  saturation  occurs  when  the  jet  deflects 
0.5bg  at  the  splitter;  hence,  a  jet  angle  of  1/40  radians  or  1.4  deg 
causes  a  pressure  differential  of  twice  the  pressure  recovery.  Typical 
null  offsets  are  2  percent.  Hence,  the  jet  deflection  is  of  the  order 
of  0.03  deg  or  0.5  mr.  Such  large  null  offsets  in  a  high  gain  system 
can  easily  cause  saturation  in  one  or  two  stages. 

For  example,  if  the  gain  is  10^  and  the  output  saturation  is 
20  kPa,  then  an  input  of  2.0  x  10~5  kPa  (0.02  Pa)  causes  saturation.  A 
rate  sensor  operating  typically  at  50  Pa  recovers  at  saturation  about 
20  Pa,  but  0.02  Pa  output  saturates  the  gain  block,  since  this  pressure 
of  0.02  Pa  represents  1/1000  of  the  saturation  pressure,  a  0.1-percent 
offset  is  intolerable.  So  one  needs  an  offset  of  less  than  at  least 
0.01  percent  or  200  times  less  than  normal  offset.  This  is  a  jet 
deflection  or  a  geometric  alignment  of  less  than  0.00015  deg  or  a 
displacement  of  the  splitter  of  0.0001bg  relative  to  the  nozzle.  For 
bg  =  0.5  ran,  this  represents  the  physically  impossible  tolerance 
requirement  of  50  ym  (2  x  10“6  in.). 


1ZF.  Manion  and  T.  Drzewiecki,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  _I,  Harry  Diamond  Laboratories  (October  1974). 
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Many  manufacturers  have  tried  to  make  the  devices  symmetric. 
Standard  format  devices  can  be  manufactured  by  exacting  techniques  such 


as  fine  blanking25  and,  for  large  production  runs,  this  appears  to  be  a 
reasonable  method.  Since  alignment  of  the  nozzle  and  the  splitter  is 
crucial,  vertically  laminated  devices  as  shown  in  figure  34  would  aid 
symmetry.  Capillary  shunts  are  provided  in  one  of  the  laminates  to  give 
temperature  compensation.  This  method  of  fabrication  closely  follows 
that  of  O'Neal,26  who  proposed  the  stacked  port  amplifiers  shown  in 
figure  35.  The  improvement  in  null  offset  is  nominally  an  order  of 
magnitude,  as  shown  in  figure  36,  where  the  null  is  referenced  to  the 
relative  rate. 


Figure  35.  Vertically  laminated  fluidic  amplifiers  (after  C.  D.  O'Neal 
U.S.  Patent  3,468,331,  23  September  1969). 
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Figure  36.  Horizontal  and  vertical  laminate  construction 

for  null  offset  versus  supply  flow:  MIL-H-5606 
hydraulic  fluid,  saturation  rate  900  deg/s,  and 
threshold  less  than  0.1  deg/s. 
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More  recently,  another  approach  was  tried.  Recognizing  that  an 
output  error  exists,  we  fed  this  signal  back  through  a  resistive  circuit 
to  the  inputs  (negative  feedback)  and  used  the  jet  deflection  gain  to 
achieve  a  net  attenuation.  At  about  the  same  time,  we  tried  a  second 
approach  that  uses  the  flow  difference  generated  between  an  orifice  and 
a  capillary  to  deflect  the  jet  in  a  direction  opposite  to  its  null 
deflection  (fig.  37).  This  technique  requires  matching  resistors  to  the 
particular  null  characteristic  at  hand,  although  variable  orifices 
(valves)  and  variable  capillaries  can  be  used  to  provide  adjustments. 

In  yet  another  technique,  we  used  the  jet  deflection  augmen¬ 
tation  principle.  Equation  (43)  relates  additional  jet  deflection  to  an 
initial  displacement  in  the  presence  of  control  edges.  When  the  bias 
control  pressure  is  some  positive  value,  the  transfer  function  is  less 
than  unity.  More  simply,  equal  bias  acts  as  a  suppressor  of  jet  deflec¬ 
tion.  This  is  intuitively  obvious  when  one  recognizes  that  deflection 
gain  decreases  with  increasing  bias.  Bias  flow  straightens  the  jet  out 
between  the  control  edges  and  gives  dramatic  feedback  results. 

O'Keeffe27  considered  the  problem  of  multistage  fluidic  ele¬ 
ments.  The  offset  of  the  first  stage  dominates  in  that  it  is  multiplied 
by  the  gain  of  the  later  stages.  By  use  of  high  resistance  feedback  in 
only  the  first  stage  (to  prevent  overall  loss  of  gain),  O'Keeffe  reduced 
overall  three- stage  LPA  offsets  of  40  percent  to  less  than  2  percent. 
More  important,  the  absolute  offset  was  kept  constant  with  supply  pres¬ 
sure  as  shown  in  figure  38.  As  a  side  benefit,  the  range  of  constant 
gain  was  improved  by  a  factor  of  two.  In  figure  39,  the  original  gain 
change  is  over  30  percent  for  a  Reynolds  number  change  of  a  factor  of 
two  (140:70),  as  is  expected  by  application  of  the  first  rule  of  thumb 
in  section  3.2.  With  null  compensation,  the  gain  change  is  less  than  15 
percent  over  the  same  100-percent  Reynolds  number  range.  The  useful 
range  of  a  10- percent  change  in  gain  is  extended  for  the  gain  block  from 
Np  of  110  to  140  to  of  85  to  140.  This  latter  is  the  same  change 

in  the  Reynolds  number  that  one  would  expect  of  an  uncompensated  nozzle 
in  pneumatic  operation  over  a  temperature  range  of  -40  to  +80  C.  Per¬ 
formance  is  considerably  superior  at  higher  aspect  ratios,  that  is,  a 
>  1.  The  amount  of  null  attenuation  by  error  feedback  is  minimal.  The 
basic  effects  are  those  of  the  bias  level  and  the  induced  differential 
flow  through  unequal  feedback  resistances.  (See  fig.  39,  p.  71.) 

The  important  contribution  of  this  research  is  that  there 
exist  a  temperature  and  a  null  compensation  with  a  single  circuit  that 
do  not  consume  any  additional  flow  or  power  and  that  do  not  depend  on 
power  supply. 

27D.  O'Keeffe,  Null  Offset  Compensation  of  Fluidic  Amplifiers  Using 
External  Negative  Feedback,  Proceedings  of  HDL  Student  Technical  Sym¬ 
posium,  Harry  Diamond  Laboratories  (August  1978). 
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Figure  37.  Null  offset  compensation  of 

laminar  jet  angular  rate  sensor 
by  orifice-capillary  control 
bias . 

Other  techniques  that  have  been  tried  include  bending,  vanes  in 
the  supply,  movable  obstructions  in  the  jet,  and  variable  in-line  resis¬ 
tors  between  the  sensor  and  the  gain  block. 

In  conclusion,  circuit  techniques  such  as  bias  control  and 
feedback  have  the  best  chance  of  universal  application  due  to  their 
simplicity  and  performance.  The  mechanical  techniques  may  tend  to  be 
susceptible  to  vibration  environments  and  be  less  reliable  due  to  moving 
parts.  Bending  in  itself  has  shown  to  be  susceptible  to  thermal  stress¬ 
ing  due  to  gradients.  In-line  resistors  yield  unequal  response  in 
different  directions.  With  better  fabrication  techniques,  null  is  not 
an  insurmountable  problem. 
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Figure  38.  Effect  of  feedback  compensation  on 
null  offset. 


3.5  Drift  and  Drift  Compensation 

Drift  is  closely  related  to  variation  of  null  offset  (sect. 
3.4).  Drift  is  essentially  the  change  of  output  level  as  a  function  of 
time.  The  time  rate  of  change  of  output,  therefore,  is  the  critical 
parameter.  Since  output  can  change  with  other  variables,  we  write  the 
total  derivative  as 


dAP  3AP  dP  3AP  3AP  dQ  3AP  , 

_ o  _  _ o  _ s  +  _ o  dT  +  _ o  _ s  +  _ o  d^ 

dt  ”  3P  dt  3T  dt  +  3Q  dt  +  3y  dt  ' 

s  s 


where  the  offset  depends  on  supply  pressure,  flow,  temperature  (T),  and 
other  causes,  x*  All  these  effects  are  Reynolds  number  effects?  hence, 
each  term  can  be  written  in  terms  of  the  Reynolds  number  as 


dAP 

_ o 

dt 


3AP  /  3N  dP  3N_  d Q  3N„  3N„  J 

_ S  I  — 5 _ 2.  +  _*_£  +  _£iI  +  _R3x 

3N_  \3P  dt  3Q  dt  3T  dt  3x  dt 

R  \  S  S 


(88) 
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If  the  kinematic  viscosities  are  the  same,  there  is  no  difference  in  net 
offset.  The  design  rule,  therefore,  is  not  to  introduce  more  gain  than 
necessary  to  overcome  the  difference  in  absolute  viscosity. 

Typically,  an  LJARS  may  exhibit  an  offset  of  2  percent  at 
100  percent  of  operating  pressure  or  3APq/3Ps  ~  0.02.  If  a  typical 
regulator  can  hold  only  ±2  percent  of  a  set  pressure,  then  3Pg/3t  * 
0.02/t.  This  value  of  3APQ/3Pg  is  related  only  to  the  inverse  of  the 
CMRR  where 


CMRR 


Change  of  Pc  at  a  constant  supply  has  the  same  effect  as 
holding  the  control  constant  and  changing  the  supply  or 


3AP  -3AP 

_ o  __  _ o 

3P  '  3  P 

c  s 


since  it  is  merely  a  change  of  bias. 

Therefore,  there  are  two  basic  means  of  controlling  drift:  re¬ 
duction  of  offset  and  reduction  of  supply  variations.  The  reduction  of 
offset  has  already  been  considered.  The  reduction  of  supply  variations 
is  not  within  the  scope  of  this  project;  however,  typical  solutions  use 
multiple  regulators  in  series,  bottled  gas,  precisely  controlled  elec¬ 
tric  motor  driven  pumps,  and  feedback  controllers. 

3. 6  Dynamic  Range  and  Maximum  Rate 

The  dynamic  range  of  a  device  is  defined  by  the  ratio  of  the 
input  signal  that  causes  saturation  or  maximum  output  to  the  input 
signal  that  causes  the  first  perceptible  output  signal.  Btor  a  device 
with  a  constant  gain,  this  is  the  same  as  the  ratio  of  saturation  output 
level  to  output  noise  level .  Ever  since  the  inception  of  laminar  flu¬ 
idic  devices,  the  concept  of  noise  has  been  of  concern.  Theoretically, 
a  laminar  flow  has  no  noise  in  and  of  itself;  hence,  the  theoretical 
dynamic  range  limit  for  an  LJARS  or  an  LPA  is  infinite.  Laboratory 
experiments  during  this  project  have  shown  that  an  LPA  can  have  dynamic 
ranges  in  excess  of  10^.  However,  to  meaningfully  establish  the  dynamic 
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range  of  an  LJARS,  one  must  be  able  to  predict  Its  maximum  rate  and 
noise. 


The  maximum  rate  that  can  be  sensed  by  an  LJARS  is  that  which 
produces  a  jet  deflection  at  the  splitter  that  gives  a  maxinum  recovered 
pressure.  Drzewiecki16  has  shown  that  the  maximum  output  of  one  port  is 
twice  the  jet  centered  pressure  recovery;  hence,  the  maximum  differ¬ 
ential  pressure  in  one  direction  is  two  times  the  jet  centered  pressure 
recovery  for  one  output  port. 


AP  =  2P 

o  max  rec 


(95) 


The  rate  sensitivity  is  defined  as 


9 

Hence,  the  rate  at  which  maximixn  recovery  is  attained  is 

2P 

9  =  . 

max  S 


(96) 


When  the  result  from  equation  (26)  is  substituted  into  equation  (96), 
the  result  is  ( in  degrees  per  second) 


9 

max 


57.3c  _vN 
_ d  R 

2b2ZO 

s 


(97) 


When  ZO  =  X2  /2  (an  efficient  nozzle  at  a  high  Reynolds  number),  we  use 
sp 

the  result  of  the  frequency  at  a  phase  shift  of  -90  deg  from  equa¬ 
tion  (65)  so  that 


•  =  16(57.3) 

max  X  90 


(98) 


sp 

This  result  indicates  that  the  maximxm  rate  is  directly  related  to  the 
frequency  response. 


16T.  M.  Drzewieckl,  Fluerics  38  s  A  Computer  Aided  Design  Analysis  for 
the  Static  and  Dynamic  Port  Characteristics  of  Laminar  Proportional 
Amplifiers,  Harry  Diamond  Laboratories  HDL-TR-1758  (June  1976). 


For  the  LJARS  considered  in  figure  12,  operating  in  hydraulic 
fluid,  an  operating  rate  sensitivity  is  2  x  10~3  psi/(deg/s).  That  same 
device  exhibits  a  bandwidth  fgg  ~  12  Hz  as  shown  in  figure  19  and  com¬ 
puted  from  equation  (65)  when  X  =  20,  bg  =  0.9  mm,  *  400,  c^  =  0.5, 
and  v  *  1.8  x  10~5  s.  Therefore,  the  maximun  rate  that  can  be 

observed  is  550  deg/s  if  the  device  is  linear  up  to  saturation.  Figure 
40  shows  the  output  signal  and  corresponding  noise  levels  for  the  same 
LJARS  in  response  to  various  levels  of  rate  input.  Hie  output  is 
significantly  identifiable  when  the  input  rate  is  as  low  as  0.1  deg/s. 
Hie  noise  level  is  roughly  one-third  that  output?  hence,  one  should  be 
able  to  discriminate  a  signal  of  one- sixth  the  input  or  about 
0.02  deg/s.  This  signal  results  in  a  dynamic  range  of  27,500  and  is  of 
the  same  order  of  the  dynamic  range  of  14,000  reported  by  Young.11 

Most  noise  cancels  out  in  the  differential  mode  by  virtue  of 
the  CMRR  of  a  device.  If  the  outputs  are  close  together,  then  the  noise 
of  the  system  is  imposed  simultaneously  on  both  and  causes  a  zero  dif¬ 
ferential  signal.  However,  there  exists  a  finite  CMRR  and  also  a  non¬ 
zero  gradient  of  output  differential  pressure  with  supply  pressure. 
Typical  values  of  9APC/3PS  =  0.02  are  discussed  in  section  3.5.  Noise 
can  be  at  any  point  of  the  frequency  spectrum;  hence,  drift  may  be 
merely  the  low-frequency  component  due  to  supply  noise.  It  follows  that 
the  minimum  signal  that  can  be  discerned  is  half  the  noise  level  divided 
by  the  sensitivity  (or  the  equivalent  rate  of  half  the  noise). 

0.  SAP 

•  o 

0  *  . . 

min  S 


noise 


(99) 


The  noise  level  due  to  a  fluctuating  supply  is  approximately 

3AP 

«  r — -  AP  , 
noise  3P  s 
s 

where  APg  is  the  amplitude  of  the  change  of  the  supply  or  the  degree  to 
which  one  can  regulate  the  supply  pressure  or  flow.  Therefore,  the 
dynamic  range,  DR,  is 


A 

max 


16(57. 3)fgoS 


Hr.  Young,  Development  of  a  laminar  Angular  Fate  Sensor  (LARS),  Pro¬ 
ceedings  of  HDL  Fluidic  State-of-the-Art  Symposium,  II,  Harry  Diamond 
laboratories  (October  1974), 
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X  — —  AP 
sp  3P  s 
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But  the  CMRR  is  merely  the  pressure  gain  divided  by  the  supply 
sensitivity  (sect.  3.5);  hence, 

1835fgQS(CMRR) 

°R  =  X  G  AP 
sp  p  s 


■■••MW 


(101) 


For  the  above  problem,  fgQ  =  12  Hz,  S  =  2  x  10~3,  X  =  20,  CMRR  ~ 
400, 23  G  ~  10,  APS  *»  0.005  psi,  and  DR  =  1.75  x  104.  The  supply 
pressure  is  about  5  psi;  hence,  the  regulation  requirement  is 
±0.1  percent  to  maintain  dynamic  range.  The  amount  of  noise  suppression 
or  supply  regulation  required  for  a  given  dynamic  range  can  be 
determined  from  equation  (101)  so  that 


1835fg()S(CMRR) 

X  G  DP 
sp  p  R  s 


/p  \ 
4  I  reel 

'°»\V 


(102) 


For  example,  to  obtain  a  dynamic  range  of  106  as  in  the  LPA  experiments, 
three  regulators  and  a  very  large  tank  are  used  in  series  to  supply  the 
device.  From  equation  (102),  for  a  typical  LJARS  to  have  the  same 
Dd  =  10^,  AP-/P-  =  1.5  x  10"5  or  a  regulation  to  ±0.0015  percent  of 
nominal  supply  pressure  over  the  entire  bandwidth  of  amplitude  response 
or  roughly  120  Hz  in  this  case.  A  tank  or  an  accumulator  can  be  tuned 
as  a  low  pass  filter  to  a  low  enough  frequency  to  take  out  all  high  fre¬ 
quencies.  Tne  regulators  act  to  take  out  low-frequency  components  if 
they  are  good  enough. 

Most  mechanical  regulators  have  a  typical  suppression  ratio  of 
about  1000.  For  example,  a  100-kPa  change  in  upstream  pressure  results 
in  a  change  of  only  0.1  kPa  at  the  load.  Hence,  two  in  series  cause  a 
suppression  of  106.  Thus,  changes  at  the  compressor,  the  pump,  or 
another  supply  can  virtually  be  eliminated.  Typically,  when  the  load 

23M.  Cycon,  T.  Sutton,  and  H.  Gamble,  Environmental  Test  Program  for 
the  Evaluation  of  Fluidic  Elements,  AiResearch  Manufacturing  Co.  of 
Arizona,  Phoenix,  AZ,  HDL-CR-7 5-187-1  (12  December  1975). 
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resistance  changes,  as  when  the  temperature  changes,  pneumatic 
regulators  do  not  respond  until  a  ±0.2-imn-Hg  change  occurs  at  best. 
Consequently,  the  regulator  appears  as  a  resistance  in  series  with  a 
variable  resistance.  If  the  load  is  predominantly  an  orifice,  as  is  the 
regulator,  then  there  is  little  relative  change  in  resistance  ratio  and 
the  supply  pressure  is  unaffected.  If,  on  the  other  hand,  the  load  is 
predominantly  a  viscous  restriction,  then  there  can  be  appreciable 
variation  in  the  pressure  division  and  the  supply  pressure  to  the 
load.  Typically,  however,  most  LPA's  and  LJARS's  may  be  considered  as 
orifices;  hence,  changes  with  temperature  should  not  change  the  supply 
pressure  significantly.  One  can  see  the  difficulty  of  using  long-nozzle 
devices . 

In  most  mechanical  regulators,  an  absolute  value  of  pressure  is 
required  to  overcome  friction,  which  is  about  ±0.2  mm  Hg  regardless  of 
the  pressure  regulated.  The  higher  the  nominal  pressure,  the  better  the 
percentage  of  regulation.  The  pneumatic  rate  sensor  and  the  amplifier 
described  in  figures  30  and  31  require  a  minimum  supply  pressure  of 
0.5  mm  Hg  to  the  rate  sensor  and  a  maximum  of  20  mm  Hg  to  the  last 
amplifier  stage.  Ttoe  pressure  regulation  can  be  kept  to  ±0.2  mm  Hg  so 
that  the  rate  sensor  supply  varies  at  most  0. 005  mm  Hg  when  the  pressure 
is  divided  through  a  series  of  resistors.  The  dynamic  range  from  equa¬ 
tion  (101)  therefore  is  greater  them  1.75  x  lo2  when  CMRR  =  400  and  the 
regulation  is  1  percent.  Fbr  a  regulation  to  ±0.1  percent  and  an 
increase  in  the  CMRR  to  4000,  the  dynamic  range  exceeds  1.75  x  105. 
Such  regulation  and  CMRR  are  not  too  difficult  to  achieve  with  present 
techniques . 

To  improve  the  response  or  the  sensitivity  by  elimination  of 
friction  due  to  moving  parts  of  a  regulator,  one  may  use  a  fluidic 
regulator.  This  is  a  high-gain  amplifier  with  negative  feedback.  The 
regulation  is  limited  by  the  available  forward  gain.  An  operational 
amplifier  saturated  to  one  side  may  be  used  to  deliver  constant  pres¬ 
sure.  The  penalty  for  such  operation  is  flow  consumption.  In  an 
air-driven  system,  flow  may  not  be  a  problem  whereas,  in  a  hydraulic 
system,  flow  is  generally  at  a  premium. 

3. 7  Power  Supply  Considerations 


We  have  already  discussed  power  supply  requirements  in  the 
sections  dealing  with  environmental  effects,  offset,  drift,  and  dynamic 
range.  It  remains  only  to  discuss  contamination  and  the  implementation 
of  the  required  power  supplies.  Unfortunately,  the  design  of  fluidic 
power  supplies  has  been  generally  neglected.  Fluidic  circuits  have  been 
expected  to  operate  and  have  operated  in  near-impossible  conditions  such 
as,  for  example,  the  DC-10  aircraft  thrust  reverser  control  that  oper¬ 
ates  off  compressor  bleed  air  in  the  jet  engine  environment. 28  However, 

28T.  Sutton  and  W.  Anderson,  Aerospace  Fluidics  Applications  and 
Circuit  Manufacture,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  V_,  Harry  Diamond  Laboratories  (October  1974), 


when  one  must  operate  a  precision  control  circuit,  as  with  a  platform 
stabilization  system,  attention  must  be  paid  to  the  power  supply.  Con¬ 
siderable  attention  has  been  paid  to  contamination  effects  and  means  of 
control .  29“32  The  basic  conclusion  that  can  be  drawn  from  these  reports 
is  that  a  laminar  fluidic  circuit  operating  on  well-filtered  air  has  a 
mean  time  to  failure  (MTTF)  of  about  25,000  hr;  a  hydraulic  circuit  has 
about  50,000  hr.  'ftiis  is  equivalent  to  over  3  years  of  continuous  duty 
and  30  years  on  a  10-percent  duty  cycle.  Failures  are  nearly  all 
attributable  to  clogging33”35  an(jf  hence,  are  reversible. 

In  general,  power  supplies  are  categorized  as  mechanical , 
chemical,  «vd  storage. Mechanical  devices  are  all  manners  of  pumps, 
compressors  fans,  and  engine  bleed.  These  are  the  typical  sources  of 
fluidic  power  for  continuous  duty  fluidic  systems.  Chemical  supplies 
require  a  reaction  of  sorts ,  and  the  gaseous  products  of  combustion  or 
reaction  are  used  as  the  working  fluid.  These  are  generally  one-shot 
supplies.  Storage  supplies  are  in  effect  fluidic  rechargeable  batter¬ 
ies.  These  include  bottled  gas,  liquefied  gas,  and  pneumatically 
pressurized  hydraulic  accumulators.  These  storage  power  supplies 
provide  the  most  uniform  quiet  flow  and  are  well  suited  to  precision 
work.  Used  with  a  mechanical  device  to  provide  intermittent  recharging, 
these  supplies  offer  the  best  solution  to  long-time  power  supply  for 
sophisticated  circuitry,  provided  that  the  flow  consumption  is  not 
excessive. 
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General  rules  of  thumb  may  be  applied  to  compute  the  flow 
consumption  of  laminar  devices.  A  pneumatic  device  with  a  nozzle  width 
and  a  throat  length  of  0.5  urn  operating  at  a  nominal  Reynolds  number  of 
120  consumes  0.3  liter/min  of  air.  Flow  at  constant  N^  is  propor¬ 
tional  to  nozzle  width,  independent  of  depth,  and  flow  is  proportional 
to  N^;  hence. 


b  N^ 

5  (b3,  hr)  -  0.3  5*5  . 

8  (V  “hW  ■  0 "i  (*=?) 


(103) 


With  this  formula  and  the  number  and  the  size  of  the  elements,  one  may 
readily  compute  the  capacity  requirements  of  any  pneumatic  laminar 
fluidic  system.  For  hydraulic  operation,  the  same  device  consumes 
6  x  10~6  m3/s;  hence, 

2  (V  NR)oil  =  10"7bs  (nnn>  %(£)  *  (104) 

In  English  units,  an  element  consumes  0.1  gpm  so  that 


®(bs'  Nr)  =  °*0017bs  (mm)  Nr  (ot®)  • 

For  comparison,  a  typical  vortex  rate  sensor  in  a  tank  gun 
stabilization  system  consumes  10  times  that  of  a  comparable  LJARS.4 

3.8  Noise  Effects 


Noise  can  arise  from  two  main  places:  the  power  supply  or  the 
vents.  From  either  source,  common-mode  pressure  noise  changes  the  oper¬ 
ating  point.  Directed  noise  at  the  vents  of  an  LJARS  can  cause  a  stand¬ 
ing  pressure  that  can  cause  jet  deflection  and  give  rise  to  a  false  or 
offset  signal.  Directed  oscillatory  noise  may  be  rectified  in  the  chan¬ 
nels  and  geometry  of  the  vent  cavity  giving  rise  to  a  dc  pressure. 
Directed  or  preferential  noise  can  occur  due  to  unequal  reflections  from 
asymmetrically  placed  reflecting  surfaces  in  pneumatic  fluidic  circuits 
that  vent  directly  into  ambient. 


UK.  Chon  and  V .  Neradka ,  Development  of  Critical  Components ,  TriTec, 
Inc.,  Columbia,  MD,  DAAG-39-77-C-0173  (July  1978). 
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No  such  noise  problems  have  been  encountered  in  liquid  oper¬ 
ation  since  the  vent  regions  are  closely  coupled  and  collected.  Close 
coupling  and  vent  collecting  in  pneumatic  circuits  also  provide  satis¬ 
factory  operation,  but  often  produce  other  problems  such  as  increased 
circuit  and  manifolding  complexity. 

Sonic  isolation19'22  is  one  method  that  has  proven  success¬ 
ful.  In  this  method,  a  separate  box  is  provided  within  which  the 
fluidic  package  is  encased.  This  box  has  a  sonic  orifice  bleed  to  the 
outside.  The  fluidic  devices  are  pressurized  to  give  the  desired  pres¬ 
sure  drop  across  them  (p  -  Pv]  ,  and  the  sonic  orifice  is  adjusted  to 
provide  a  two-to-one  pressure  ratio  from  the  box  to  ambient,  insuring 
choked  conditions.  This,  incidentally,  causes  the  density  of  the  air  to 
double.  If  the  box  is  acoustically  insulated,  the  system  should  be 
adequately  isolated.  More  typically,  external  disturbances  such  as 
movement  of  solid  objects  may  also  cause  dramatic  output  changes  as 
shown  by  Young9?  and  reproduced  in  figure  41.  The  effects  of  motion  can 
be  eliminated  satisfactorily  when  a  foam  material  covers  the  vents.  The 
foam  also  acts  as  an  acoustic  baffle  and  reduces  noise  input.  High- 
frequency  hissing  is  a  typical  noise  source  and  is  adequately  eliminated 
by  foam  rubber  surrounding  the  vents.  Care  should  be  exercised  so  that 
the  vent  is  not  adversely  affected. 

Noise  or,  more  specifically,  selected  frequency  acoustic 
signals  can  cause  a  laminar  jet  that  is  near  transition  to  go  into 
transition.  Specifically,  a  pneumatic  WARS  operating  at  oNR  ~  1300  was 
found  to  go  into  transition  and  develop  an  offset  when  subjected  to 
acoustic  signals  of  about  100  dB  at  between  150  and  250  Hz  when  a  nearby 
pump  was  turned  on.  This  action  strongly  indicates  that  one  must  be 
sure  to  operate  at  cjNr  <  1000  as  stipulated  in  section  3.1.  In  the 
above  example,  when  the  supply  pressure  was  reduced  20  percent,  the 
effect  of  the  pump  noise  was  eliminated.  Its  elimination  suggests 
that  oNr  -  1000  is  about  the  limit  of  absolute  laminar  jet  stability. 
As  the  Reynolds  number  was  decreased,  as  expected,  higher  frequencies 
were  needed  to  set  off  transition  until  it  did  not  occur  at  all.  Its 
absence  suggests  a  research  area  on  determination  of  the  stability  of 
confined  laminar  jets.  Regions  of  stability  with  frequency  (wave 
number),  amplitude  of  perturbation,  and  Reynolds  number  should  be 
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determined  in  something  akin  to  neutral  stability  curves  as  determined 
for  laminar  boundary  layers. 


Figure  41.  Oscilloscope  traces  of  pneumatic  laminar  jet  angular  rate 
sensor  output  with  uncovered  vents:  (a)  still  air  and 
(b)  ambient  air  disturbed  by  waving  of  hand; 
output  =  0.025  mm  Hg/cm,  Pace  pressure  transducer  =  0  to 
0.03  psid,  and  supply  pressure  =  12  mm  Hg  (no  rotation) . 

3.9  Operating  Fluids 

It  has  been  tacitly  assumed  that  the  fluids  used  are  New¬ 
tonian.  Non-Newtonian  fluids  may  be  used;  however,  device  operation 
will  be  somewhat  different.  A  fluid  that  shears  readily  will  cause 
higher  pressure  recovery;  hence,  the  sensitivity  may  be  greater  since 
frictional  losses  can  be  greatly  reduced.  On  the  other  hand,  temper¬ 
ature  characteristics  of  such  liquids  generally  are  not  well  behaved,  so 
it  is  probably  wise  to  stay  clear  of  these  and  thixotropic  materials. 

In  general,  there  is  no  fundamental  difference  in  the  operation 
of  an  LJARS  when  liquids  or  gases  are  used.  This  paper  considers  the 
most  common  fluids,  oil  and  air.  Many  fluidic  circuits  are  operated 
with  water  as  well,  but  water  is  not  advisable  for  operation  over  the 
military  temperature  range.  However,  in  aqueous  environments  such  as  a 
torpedo  guidance  system,  the  local  available  fluid,  sea  water,  would 
pose  no- theoretical  limitations. 

The  big  difference  in  operating  with  liquid  working  fluids  is 
the  additional  manifolding  required  to  return  the  flow  to  sump  if 
necessary.  Extreme  care  must  be  exercised  to  closely  cross  couple  the 
vent  regions  on  both  sides  of  the  jet  to  prevent  pressure  gradients  from 
forming  and  deflecting  the  jet.  Early  attempts  to  make  a  hydraulic 
LJARS  resulted  in  unacceptable  offsets  that  were  later  relieved  when 


cross  coupling  was  provided.  Figure  42  shows  a  typical  lamination  of  an 
LJARS  used  in  this  program.  Figures  43  and  44  show  the  plates  required 
to  provide  cross  coupling.  The  jet  bounding  plate  forms  the  bounding 
plane  for  the  jet,  and  the  vent  cavity  plate  couples  the  vent  regions  on 
both  sides  of  the  jet. 


NOTE:  RATE  SENSOR  DIMENSIONS  ARE  SUPPLY  WIDTH  (B,);  B,  =  0.020. 


Figure  42.  Typical  laminar  jet  angular  rate  sensor 
horizontal  lamination. 

The  cavity  provides  a  low-impedance,  uniform  pressure  ground 
that  cannot  support  lateral  pressure  gradients.  This  is  typically 
called  a  floating  ground  since  the  vent  pressure  level  depends  on  the 
impedance  back  to  sump  and  the  flow  through  the  device.  It  is  usually 
advantageous  to  provide  this  type  of  venting  on  both  the  top  and  the 
bottom  of  the  LJARS  itself.  In  this  manner,  vertical  gradients  can  be 
minimized;  however,  coupling  on  one  side  only  may  be  adequate  for  most 
applications.  When  vertical  laminations  are  used,  the  cross  coupling  is 
also  a  practical  necessity.  Figure  34  shows  a  vertically  laminated 
device  with  oval  cutouts  on  top  and  bottom.  These  form  the  cross- 
-coupling  cavity,  and  the  outer  laminations  provide  the  vent  commun¬ 
ication  to  these  cavities.  The  flow  returning  to  sump  must  be  collected 
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from  a  point  of  symmetry  so  that  an  asymmetrical  flow  distribution  is 
not  formed.  The  flow  patterns  (fig.  45)  may  indeed  deflect  the  jet, 
especially  since  we  are  expecting  to  distinguish  a  jet  deflection  on  the 
order  of  10“6  deg.  In  the  asymmetric  case  (fig.  45),  the  flow  to  the 
lower  vent  area  must  be  turned  upward;  hence,  a  pressure  develops  that 
deflects  the  jet  to  give  a  spurious  signal.  When  vertical  laminates  are 
used,  multiple  side  outlets  must  be  equally  spaced. 


vent  areas. 


Figure  44.  Cross-coupling  vent  cavity  plate. 


Figure  45.  Flow  patterns  with  collected  vents. 


3. 10  Fabrication 

Most  fluidic  devices  are  made  of  laminations  to  provide  for  a 
cancellation  effect  of  mechanical  errors  when  they  are  stacked. 
Sutton28  discusses  typical  industrial  methods  of  photochemical  etch¬ 
ing.  Phillippi25  discusses  fine  blanking,  which  is  a  high-precision 
stamping  process.  It  has  the  great  advantage  of  superior  repeatability 
from  laminate  to  laminate  once  a  die  has  been  custom-tailored  to  give 
the  desired  geometry.  Custom-tailoring  leads,  however,  to  high  die 
cost.  To  amortize  the  die  cost  so  that  laminates  are  competitively 
priced  with  those  of  metal  etching,  many  thousands  of  pieces  must  be 
made.  Interestingly,  laminations  can  be  thickened  so  that  the  cost  per 
unit  thickness  can  guickly  be  competitive. 


25P.  M.  Phillippi ,  A  Study  of  Fineblanking  for  the  Manufacture  of 
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Elements  produced  by  electroforming  have  the  disadvantage  of  a 
high  cost  die  and  a  geometrical  limitation  on  element  depth.  Normally, 
devices  are  limited  to  an  aspect  ratio  of  1. 

Elements  may  be  fabricated  in  plastic  by  molding  or  casting. 
Normally,  military  fluidic  hardware  is  not  plastic  due  to  lack  of 
ruggedness  and  dimensional  stability.  Dimensional  accuracy  also  is 
more  difficult  to  achieve  since  the  three-dimensional  geometry  is  made 
at  once. 


A  relatively  new  fabrication  technique  is  electrical  discharge 
machining  (EDM).  A  wire  cuts  out  the  desired  shape  in  a  thick  piece  of 
material.  Tolerance  can  be  maintained  by  computer  control.  In  this 
manner,  single  laminations  of  the  desired  aspect  ratio  can  be  manu¬ 
factured.  Hie  cost  per  laminate  is  prohibitively  high,  up  to 
$25/laminate  or  very  nearly  100  times  more  them  with  etching.  This  cost 
is  essentially  equivalent  to  that  of  optical-tracing  milling  machined 
elements  since  the  same  computer  control  is  required.  These  techniques 
are  used  for  small  orders  and  laboratory  tests. 

Typical  costs  of  production  items  are  $0. 25/laminate  for 
metal-etched  laminations.  For  a  device  of  5  1  2.6,  bg  =  0.5  mm,  and  a 
laminate  thickness  of  0.1  mm  for  the  LJARS  and  0.25  mm  for  auxiliary 
manifolds,  13  LJARS  laminates  are  required  to  make  o  =  2.6.  Two 
bounding  plates  (fig.  43)  are  needed,  and  the  cavity  depth  is  twice  the 
height  of  a  device  as  a  rule  or  2.6  mm  in  this  case.  This  depth 
requires  10  cavity  laminations  on  each  side.  To  complete  the  device  for 
circuit  integration,  a  cover  and  a  baseplate  are  required.  Thirty-seven 
laminates  cost  $9.25  at  production  prices  not  including  labor  and 
bonding  charges.  At  research  prices  of  up  to  $2  each,  37  laminates  cost 
up  to  $74.  This  is  well  below  $2000  or  more  than  a  gyroscope  of  similar 
functional  properties  may  cost. 

3.11  Testing  Procedures 

The  verification  of  the  operational  characteristics  of  any 
device  is  of  paramount  importance  in  any  practical  application.  For 
this  reason,  we  look  at  the  methods  required  for  obtaining  valid  experi¬ 
mental  data.  The  characteristics  of  an  LJARS  are  not  unlike  those  of  an 
LPA.  Indeed,  the  pressure  gain  of  the  device  is  one  of  the  quantities 
required  to  determine  dynamic  range  as  shown  in  equation  (101);  hence, 
the  tests  must  include  all  those  required  of  an  LPA  with  the  addition  of 
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a  rate  sensitivity  test. 19,23,38  This  section  briefly  describes  the 
testing  procedures  typically  conducted  on  an  UARS  • 

Eight  critical  tests  are  required  of  an  LJARS: 

Supply  characteristic. — The  supply  flow  is  measured  against 
pressure  with  all  ports  open  and  then  blocked.  Normally  if  the  device 
is  well  vented,  there  is  little  difference  in  whether  the  ports  are  open 
or  blocked.  The  supply  pressure  is  the  pressure  drop  from  the  plenum  to 
the  vent  cavity  (or  to  ambient  if  so  vented)  .  A  supply  Reynolds  number 
is  defined  by  this  drop,  that  is, 


Pressure  recovery  characteristic. — Blocked  output  pressure  is 
measured  versus  supply  pressure.  Five  curves  are  typically  taken:  one 
at  each  of  two  outputs  with  the  opposite  output  blocked,  one  at  each  of 
two  outputs  with  the  opposite  output  open,  and  one  with  both  ports  fed 
together  to  given  an  average  value.  This  last  is  generally  used  for  the 
pressure  recovery.  The  first  four  curves  give  measures  of  the  offset 
and  the  load  sensitivity.  For  example,  both  blocked  should  not  affect 
jet  position.  One  open  may  cause  jet  deflection  and  gives  a  measure  of 
loading  effects.  This  test  may  be  repeated  for  open  or  blocked 
controls . 

Null  offset  characteristic. — Differential  blocked  output 
pressure  is  measured  versus  supply  pressure.  Curves  are  taken  with  open 
controls,  with  blocked  controls,  and  with  any  desired  compensation 
circuitry.  These  curves  give  APQ/Pg,  which  is  required  for  dynamic 
range,  drift,  environmental  effects,  and  CMRR  computations. 

Input  resistance  characteristic. — Control  flow  is  measured 
versus  control  pressure,  where  control  pressure  is  referenced  to  the 
vent  cavity.  Two  sets  of  curves  are  taken.  The  first  is  called  jet 
centered  resistance,  which  is  a  measure  of  the  control  channel  plus  the 


19G. Roe,  Fluidic  Laminar  Angular  Rate  Sensor  Research  Program, 
McDonnell  Douglas,  Titusville,  FL,  HDL-CR-7 5-020-1  (November  1975). 

23M.  Cgcon,  T.  Sutton,  and  H.  Gamble,  Environmental  Test  Program  for 
the  Evaluation  of  Fluidic  Elements,  AiResearch  Manufacturing  Co.  of 
Arizona,  Phoenix,  AZ,  HDL-CR-75-187-1  (12  December  1975). 

38Fluerics  Testing  Procedures,  Department  of  Defense  MIL-STD-1361 
(1971). 
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jet  edge  resistance  (rc  +  Ry).  This  is  obtained  by  measuring  the  total 
flow  through  a  tee  to  both  controls  simultaneously  versus  the  port 
pressure.  A  plot  is  made  of  half  that  flow  versus  pressure.  The 
procedure  is  repeated  for  the  jet  turned  off.  This  curve  gives  the 
control  channel  resistance,  Rc,  and  the  difference  gives  the  vent 
resistance,  Rv,  required  in  augmented  jet  deflection  computations.  The 
second  set  of  curves  is  control  flow  versus  pressure  with  the  opposite 
side  control  open  to  ambient.  This  with  the  jet  centered  curve  can  be 
used  to  get  the  deflection  resistance  if  necessary. 

Output  resistance  characteristic . — Output  flow  is  measured 
versus  output  pressure  (referenced  to  the  vent  cavity).  This  is  tested 
in  the  same  manner  as  the  jet  centered  input  characteristic.  The 
outputs  are  fed  together,  and  the  pressure  is  varied  from  blocked 
pressure  to  zero  (or  less  if  one  desires  to  see  vhat  happens  with  a 
suction  load)  while  the  flow  is  monitored.  A  plot  is  made  of  half  this 
flow  (the  average  of  one  output)  versus  output  pressure.  This  test  may 
be  repeated  for  the  individual  output  with  opposite  outputs  open  and 
then  blocked  to  give  a  measure  of  load  dependence.  The  inverse  slope 
of  the  output  curve  at  the  operating  point  gives  the  operating  output 
resistance,  1^,  required  to  compute  the  gain  (sensitivity)  degradation 
when  loading  into  a  circuit.  The  operating  point  gives  the  value  of 
bias  into  the  circuit.  The  gain  is  degraded  in  a  circuit  by  the  factor 
Rload/(Ro  +  Rload) '  where  Rload  is  the  inPut  resistance  of  the  circuit 
into  which  the  LIARS  feeds.  Hie  point  at  which  the  flow  is  zero  gives 
the  pressure  recovery;  the  point  at  which  the  pressure  is  zero  gives  the 
flow  recovery. 

Pressure  gain  or  pressure  transfer  characteristic . — Blocked 
output  differential  pressure  is  measured  versus  differential  input 
pressure  with  supply  pressure  held  constant.  The  slope  of  this  curve  at 
zero  differential  pressure  gives  the  pressure  gain,  G  ,  and  the  differ¬ 
ential  output  gives  the  offset.  The  linear  portion  of  this  curve 
defines  the  linear  range  of  operation  of  the  LJARS.  Also,  Gp  is 
required  for  computation  of  dynamic  range  and  CMRR. 

Sensitivity  or  rate  transfer  characteristic. — Blocked  output 
differential  pressure  is  measured  versus  applied  angular  rate  with 
supply  pressure  held  constant.  This  is  obtained  by  use  of  a  rate  table 
to  provide  input.  The  slope  of  this  curve  is  sensitivity. 

Frequency  response  characteristic .—Ratio  of  amplitudes  of 
blocked  output  differential  pressure  to  an  input  rate  signal  is  measured 
versus  frequency.  Typically,  this  is  called  a  Bode  plot  when  the  phase 
shift  is  included.  The  phase  shift  portion  is  critical  because  it  is 
used  to  define  the  bandwidth  of  the  LJARS.  A  gain-phase  meter  and  a 
servo-controlled  rate  table  are  generally  required  to  make  this 
measurement  with  minimum  effort. 


The  last  five  tests  are  conducted  at  various  supply  pressures 
or  Reynolds  numbers  within  the  operating  range  so  that  auxiliary  curves 

such  as  sensitivity  versus  Reynolds  number  may  be  obtained.  Figures  46 

to  51  show  the  experimental  data  for  the  LJARS  shown  in  figures  42  to  44 
operating  with  air,  bg  =  0.75  mm,  and  a  =  1.67.  Not  all  the  curves 
specified  above  are  presented  in  these  figures;  however,  all  pertinent 
data  are  presented.  The  Bode  plot  is  shown  in  figure  19  for  oil.  Since 
the  kinematic  viscosities  of  air  and  oil  are  almost  the  same  at  room 
temperature,  the  frequency  response  is  similar.  In  figure  51,  we  have 
replotted  the  data  of  figure  11,  for  completeness.  Figure  50  is  the 

plot  of  the  actual  data  from  which  figure  51  was  plotted.  The  uncom¬ 

pensated  outputs  in  figure  50  have  a  null  offset  on  the  order  of  2.0 
Pa.  At  a  sensitivity  of  0.05  Pa/( deg/s),  this  corresponds  to  an  error 
rate  of  40  deg/s,  so  the  LJARS  must  be  compensated.  Figure  52  shows  the 
null  offset  changed  by  the  negative  feedback  through  valves  from  output 
to  input.  There  exists  a  fairly  wide  range  over  which  the  null  can  be 
eliminated.  In  this  range  one  would  operate  the  UARS. 


Figure  46.  Experimental  supply  and  control  characteristics 
for  laminar  jet  angular  rate  sensor:  normalized 
splitter  length  =  20,  supply  width  =  0.75  mm, 
normalized  output  width  =  1.7,  and  aspect 
ratio  =  1.667. 


^—PRESSURE  RECOVERY 


Figure  47.  Experimental  pressure  recovery  and  null  offset 
for  laminar  jet  angular  rate  sensor. 


Westerman  and  Wright22  conducted  an  optimization  program  based 
on  the  single  parameter  of  sensitivity.  Figure  53  shows  a  sample  of 
silhouettes  used  in  the  program,  and  figure  54  shows  their  sensi¬ 
tivity.  Unit  FM4  corresponds  with  the  device  of  figure  42,  and  the 
sensitivity  is  a  good  order  of  magnitude  lower  than  for  FM5.  In  FM5, 
the  positive  feedback  effect  of  the  control  ports  is  used.  Interest¬ 
ingly,  the  theory  developed  in  section  2  more  than  adequately  predicts 
the  performance  of  EW4  and  the  others.  However,  all  the  devices  using 
control  edge  feedback  depend  on  the  Reynolds  number,  we  have  opted  for 
a  trade-off  with  sensitivity  for  operating  range.  If  one  desires  a 

single  high-gain  element  and  is  willing  to  pay  the  costs  of  keeping  the 
Reynolds  number  constant,  then  FW5  is  desirable. 

The  operating  characteristics  of  the  HDL  WARS  shown  in  figures 
42  to  44  are  given  in  table  1 . 


22IV.  J.  Westerman  and  R.  E.  Wright,  Evaluation  Program — Fluidic 
Laminar  Rate  Sensor — Final  Report,  McDonnell  Douglas,  Titusville,  FL, 
NAVAIR  L0252  (8  March  1974). 
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Figure  48.  Output  characteristic  of  experi¬ 
mental  laminar  jet  angular  rate 
sensor. 


Figure  49.  Experimental  pressure 
gain  of  laminar  jet 
angular  rate  sensor. 
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Figure  50.  Experimental  sensitivity  Figure  51.  Experimental  sensitivity  as 

of  laminar  jet  angular  function  of  supply  pressure, 

rate  sensor: 

(a)  blocked  controls  and 

(b)  open  controls. 
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Figure  54.  Sensitivity  versus  Reynolds  number 
of  laminar  jet  angular  rate  sensors 
of  figure  53  (from  G.  Roe,  McDonnell 
Douglas  HDL-CR-75-020-1,  November  1975) 


TABLE  1.  TECHNICAL  DATA  FOR  LAMINAR  JET  ANGULAR  RATE  SENSOR  OF  HARI  Y  DIAMOND  LABORATORIES 


supply  width, 
height . 

supply  pressure . 
aspect  ratio ,  h/b 
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4.  DISCUSSION  AND  CONCLUSIONS 


We  have  presented  an  analysis  and  a  few  rules  of  thumb  for  the 
design  of  a  fluidic  LJARS.  Throughout  this  report,  most  emphasis  has 
been  on  a  single  component  thought  to  be  representative — namely,  the  one 
depicted  in  figure  42.  The  raw  data  for  the  performance  are  given  in 
figures  46  to  51.  This  LJARS  is  not  optimal  because  optimality  is  a 
function  of  what  is  desired.  Nevertheless,  it  has  most  of  the  good 
characteristics  attributable  to  an  LJARS,  such  as  high  sensitivity,  low 
output  impedance,  ease  of  compensation  for  null  and  temperature,  and 
generally  low  sensitivity  to  temperature  changes.  Although  these 
characteristics  are  not  final  or  conclusive,  they  show  what  can  be  done 
with  minimum  effort.  The  negative  feedback  offset  compensation  yields 
satisfactory  results  that  depend  on  the  type  of  feedback  resistor.  In 
general,  the  shape  of  the  offset  curve  (fig.  52)  and  the  magnitude  of 
the  offset  may  be  affected,  thus  giving  considerable  flexibility. 

Temperature  compensation  is  more  dramatically  incorporated  in  a 
hydraulic  LJARS,  in  which  the  density  does  not  change,  than  in  a 
pneumatic  UkRS,  in  which  it  does.  In  a  hydraulic  LJARS,  sensitivity 
cam  be  decreased  with  an  increase  in  temperature,  but  in  a  pneumatic 
LJARS,  the  best  that  can  be  achieved  is  a  reduction  in  the  increase. 


The  threshold  depends  on  the  available  instrumentation  and  the 
environmental  noise.  Open- loop  measurements  indicate  thresholds  of  less 
than  0.1  deg/s  (as  fast  as  a  minute  hand  turns  on  a  clock) .  Closed-loop 
performance  on  a  pneumatic  fluidic  gun  stabilization  system  for  the 
M60/M48  tank  showed  disturbance  suppression  of  40  dB  for  signals  on  the 
order  of  1  deg/s,  which  means  that  the  rate  sensor  is  moving  at  0.01 
deg/s.  Closed-loop  operation  is  less  susceptible  to  common-mode  noise. 


To  facilitate  design,  the  important  equations  for  LJARS  operation 
are  these: 
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In  conclusion,  this  report  has  shown  how  an  LJARS  may  successfully 
be  modelled  and  has  presented  guidelines  to  build  it. 

The  applications  of  an  LJARS  fall  into  several  categories  from 
simple  rate  limit  switches  to  rate  gyroscope  replacement.  The  broadest 
applications  appear  to  be  in  guidance  and  disturbance  suppression  sys¬ 
tems  .  The  overwhelming  asset  of  the  LJARS  is  the  low  cost.  Current 
production  rate  costs  of  $0. 25/laminate  yield  cost  estimates  between  $10 
and  $100  for  an  LJARS.  A  trade-off  can  be  made  between  accuracy  (null 
offset  and  drift)  and  low  cost  when  an  LJARS  is  compared  with  a  gyro¬ 
scope.  The  cost  of  the  LJARS  increases  with  the  power  supply 

complexity. 

Another  asset  is  the  instant-on  feature.  Rie  sensor  works  as  soon 
as  a  jet  is  developed— one  jet  transport  time.  Typically,  for  a 

pneumatic  LJARS  with  a  jet  velocity  of  10  a/ a  over  a  distance  of  15  ran, 
the  LJARS  starts  in  only  1.5  ms.  This  time  can  be  compared  with  the 
several  seconds  required  to  spin  up  a  mechanical  gyroscope.  This  asset 
is  of  particular  interest  to  those  who  design  seat  ejection  systems  that 
must  be  on  line  in  less  than  10  ms  from  storage  to  provide  adequate 
performance. 

The  LJARS  appears  to  apply  to  many  guidance  and  control  systems. 
Lest  this  impression  be  imprinted  too  strongly,  one  must  remember  this 

important  caveat:  The  LJARS  is  only  as  good  as  its  power  supply. 

Reliable,  accurate,  fluidic  power  supplies  are  more  difficult  to  make 
with  the  current  state  of  the  art  than  electronic  ones,  for  which  far 
more  research  dollars  have  been  expended. 

In  addition,  the  problem  of  interfacing  must  yet  be  addressed. 
Interfacing  with  electronic  systems  involves  a  pressure  transducer, 
which  often  is  the  weak  link.  Interfacing  with  mechanical  systems 
involves  servo  actuators  or  servo  valves.  Fluidic  first- stage  servo 
valves  are  on  the  market.  Hence,  interfacing  may  not  present  a  big 
problem,  although  it  involves  moving  parts. 

We  suggest  that  the  interfacing  and  the  power  supplies  be  devel¬ 
oped.  Null  offset  compensation  and  dynamic  effects  need  intensive 
analytical  research  to  go  with  current  results  based  on  the  heuristic 
arguments  presented  here. 


LITERATURE  CITED 


1.  A.  Ostdiek,  Viscous  Vortex  Rate  Sensor,  Harry  Diamond  Laboratories 
HDL-TR-1555  (November  1971). 

2.  J.  M.  Kirshner,  A  Survey  of  Sensors — Part  II,  Proceedings  of  HDL 
Fluidic  State-of-the-Art  Symposium,  II,  Harry  Diamond  Laboratories 
(October  1974). 

3.  J.  0.  Hedeen,  Performance  and  Application  of  the  Vortex  Rate 
Sensor,  Proceedings  of  HDL  Fluidic.  State-of-the-Art  Symposium,  II, 
Harry  Diamond  Laboratories  (October  1974). 

4.  K.  Chen  and  V.  Neradka,  Development  of  Critical  Components,  TriTec, 
Inc.,  Columbia,  MD,  DAAG-39-77-C-0173  (July  1978). 

5.  H.  Ziebolz,  Characteristics  of  Hydraulic  and  Pneumatic  Relays  as 
Energy-Converting  Devices,  Instruments,  _15  (September  1942). 

6.  T.  Reader,  Laminar  Fluidic  Devices,  U.S.  Patent  3,662,772  (16  May 
1972). 

7.  A.  G.  Moore,  Angular  Movement  Sensing  Device,  U.S.  Patent  3,500,691 
(17  March  1970). 

8.  T.  M.  Rankin,  A.  G.  Moore,  and  W.  C.  Schuemann,  Fluidic  Angular 
Rate  Sensor— A  Replacement  for  Rate  Gyroscopes?  Johns  Hopkins 
Applied  Physics  Laboratory  Technical  Digest  (March-April  1969). 

9.  T.  Reader,  T.  Shaffer,  J.  Goren,  and  R.  Kantola,  Feasibility 
Investigation  of  a  Laminar  Rate  Sensor  (LARS),  Phase  I  Report, 
General  Electric  Co.  Reentry  and  Environmental  Systems,  Schen¬ 
ectady,  NY,  69SD698  (June  1969). 

10.  P.  L.  Jacobs,  Theoretical  Analysis  of  a  Two-Dimensional  Laminar  Jet 
Rate  Sensor,  Proceedings  of  HDL  Fluidic  State-of-the-Art  Symposium, 
II,  Harry  Diamond  Laboratories  (October  1974). 

11.  R.  Young,  Development  of  a  Laminar  Angular  Rate  Sensor  (LARS), 
Proceedings  of  HDL  Fluidic  State-of-the-Art  Symposium,  II,  Harry 

Diamond  Laboratories  (October  1974). 

\ 

12.  F.  Manion  and  T.  Drzewiecki,  Analytical  Design  of  Laminar  Propor¬ 
tional  Amplifiers,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  _I,  Harry  Diamond  Laboratories  (October  1974). 


98 


LITERATURE  (Cont'd) 


13.  T.  M.  Drzewiecki,  0.  N.  Wormley ,  and  P.  H.  Man. ion ,  Compute r - Ai ded 
Design  Procedure  for  Laminar  Fluidic  Systems,  J.  Dyn.  Syst.  Meas. 
Control,  97,  Series  6  (December  1975). 

14.  T.  M.  Drzewiecki,  Interpretation  of  Surface  Static  Pressure  Distri¬ 
butions  in  Fluid  Amplifier  Applications,  Harry  Diamond  Laboratories 
HDL-TR-1627  (July  1973). 

15.  H.  Schlichting,  Boundary  Layer  Theory,  McGraw-Hill  Book  Co.,  New 
York  (1960). 

16.  T.  M.  Drzewiecki,  Fluerics  38:  A  Computer  Aided  Design  Analysis 
for  the  Static  and  Dynamic  Port  Characteristics  of  Laminar  Propor¬ 
tional  Amplifiers,  Harry  Diamond  Laboratories  HDL-TR-1758  (June 
1976) . 

17.  P.  M.  Manion  and  G.  Mon,  Fluerics  33:  Design  and  staging  of  Lam¬ 
inar  Proportional  Amplifiers,  Harry  Diamond  Laboratories 
HDL-TR— 1 608  (September  1972). 

18.  J.  O.  Hinze,  Turbulence,  McGraw-Hill  Book  Co.,  New  York  (1959). 

19.  G.  Roe,  Fluidic  Laminar  Angular  Rate  Sensor  Research  Program, 
McDonnell  Douglas,  Titusville,  FL,  HDL-CR-75-020-1  (November  1975). 

20.  T.  M.  Drzewiecki,  Fluidic  Resistive  Element,  U.S.  Patent  3,980,103 
(14  September  1976). 

21.  George  Mon,  Flueric  Laminar  Gain  Blocks  and  an  operational  Ampli¬ 
fier  Scaler,  Harry  Diamond  Laboratories  HDL-TR- 1730  (December 
1975) . 

22.  w.  J.  Westerman  and  R.  E.  Wright,  Evaluation  Program — Fluidic 
Laminar  Rate  Sensor— Fined  Report,  McDonnell  Douglas,  Titusville, 
FL,  NAVAIR  L0252  (8  March  1974). 

23.  M.  Cycon,  T.  Sutton,  and  H.  Gamble,  Environmental  Test  Program  for 
the  Evaluation  of  Fluidic  Elements,  AiResearch  Manufacturing  Co.  of 
Arizona,  Phoenix,  AZ,  HDL-CR-75- 187-1  (12  December  1975). 

24.  V.  Bravo,  Feasibility  Investigation  of  a  Laminar  Rate  Sensor 
(LARS),  Phase  III  Report,  General  Electric  Co,  Schenectady,  NY, 
71SD2118  (December  1971). 

25.  R.  M.  Phillippi,  A  Study  of  Fineblanking  for  the  Manufacture  of 
Flueric  Laminar  Proportioned  Amplifiers,  Harry  Diamond  Laboratories 
HDL-TM-77-8  (May  1977). 


99 


LITERATURE  (Oont'd) 


26.  C.  D.  O'Neal,  Stacked  Port  Fluidic  Amplifier,  U.S.  Patent  3,468,331 
(23  September  1969). 

27.  D.  O'Keeffe,  Null  Offset  Compensation  of  Fluidic  Amplifiers  Using 
External  Negative  Feedback,  Proceedings  of  HDL  Student  Technical 
Symposium,  Harry  Diamond  Laboratories  (August  1978). 

28.  T.  Sutton  and  W.  Anderson,  Aerospace  Fluidics  Applications  and 
Circuit  Manufacture,  Proceedings  of  HDL  Fluidic  State-of-the-Art 
Symposium,  V,  Harry  Diamond  Laboratories  (October  1974). 

29.  W.  Westerman  and  R.  Wright,  Fluidic  Reliability  Program,  McDonnell 
Douglas,  Titusville,  FL,  L0242  (21  December  1973). 

30.  W.  Westerman,  R.  Wright,  and  G.  Roe,  Investigation  of  the  Effec¬ 
tiveness  of  Cyclone  Separators  on  Fluidic  Power  Supplies,  McDonnell 
Douglas,  Titusville,  FL,  L0243  (21  December  1973). 

31.  L.  Banaszak  and  W.  Posingies,  Hydrofluidic  Stability  Augmentation 
System  (HYSAS)  Operational  Suitability  Demonstration,  Honeywell, 
Minneapolis,  MN,  USAAMRDL-TR-77-31  (October  1977). 

32  W.  Westerman  and  R.  Wright,  Reliability  Testing  of  Laminar  Jet 
Fluidic  Elements,  McDonnell  Douglas,  Titusville,  FL,  L0266 
( 1  June  1974 ) . 

33.  R.  Comparin,  H.  Moses,  and  E.  Rowell,  Contamination  Effects  in  a 
Laminar  Proportional  Amplifier,  Virginia  Polytechnic  Institute, 
Blacksburg,  VA,  HDL-TR-74- 175-1  (June  1974). 

34.  H.  Moses  and  E.  Sancaktar,  The  Effect  of  Geometric  Changes  Due  to 
Contamination  on  Fluid  Amplifier  Performance,  Virginia  Polytechnic 
Institute,  Blacksburg,  VA,  HDL-CR-75-192-1  (November  1975) . 

35.  R.  Martin,  J.  Denneny,  and  T.  Sutton,  Three-Phase  Program  to  Define 
Requirements  for  Operation  of  Fluidics  on  Compressor  Bleed  Air, 
AiResearch  Manufacturing  Co.  of  Arizona,  Phoenix,  AZ,  76-411328 
( August  1 976 ) . 

36.  W.  westerman  and  G.  Roe,  Classification  of  Fluidic  Power  Supplies, 
McDonnell  Douglas,  Titusville,  FL,  L0244  (21  December  1973). 

37.  R.  Young,  Feasibility  Investigation  of  a  Laminar  Angular  Rate 
Sensor,  Phase  IV  Final  Report,  General  Electric  Co.  Reentry  and 
Environmental  Systems,  Schenectady,  NY,  73SD2147  (August  1973). 

38.  Fluerics  Testing  Procedures,  Department  of  Defense  MIL-STD-1361 
(1971). 


100 


NOMENCLATURE 


a  j  Net  entrainment  coefficient 

A  Cross-sectional  area  (m2) 

b  Width  dimension  (m) 

B  Normalized  width  (b/bs) 

c^  Flow  discharge  coefficient  =  Qa/Q^ 

i 

c0  Momentum  flux  discharge  coefficient  =  Ja/J^  ] 

C  Loss  coefficient 

Dr  Dynamic  range 

f  Bandwidth  (Hz) 

f ( x)  Center-line  decay  function 

G  Forward  gain 

h  Height  (m) 

H  Feedback  gain 

2 

J  Momentum  flux  =  pAU 

k  Constant 

K  Kinematic  momentum  flux  =  J/p 

Orifice  coefficient 
m  Mass  flow  rate  (kg/s) 

n  Integer 

NR  Reynolds  number  =  pbgUg/u 

P  Pressure  (Pa)  j 

,  | 

Q  Volumetric  flow  rate  (m3/s)  ' 
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NOMENCLATURE  (Oont'd) 


Radial  position  (m) 

Resistance  based  on  Q,  R  *  AP/Q  (kg/m^s) 

Resistance  based  on  m,  R1  =  AP/m  (l/m*s) 

Normalized  angular  rate  sensitivity 
Temperature  (C  or  K) 
x-direction  velocity  (m/s) 

Mean  or  reference  velocity  (m/s) 

Lateral  velocity 
Velocity  (m/s) 

Coordinate  dimension  measured  from  nozzle  exit  (m) 
Normalized  length  in  x-direction 
Lateral  deflection  distance  (m) 

Half  square  of  effective  normalized  splitter  distance 
Angle 

Boundary  layer  thickness 
Displacement  thickness 
Momentum  thickness 
Angular  rate 
Wavelength 

Absolute  viscosity  (kg/m*s) 

Kinematic  viscosity  (m2/s) 

Density  (kg/m^) 


NOMENCLATURE  (Cont'd) 


a  Aspect  ratio,  h/b_ 

T  Transit  time  (s) 

$  Phase  angle  (deg) 

X  Other  causes 

Subscripts 

a  Augmented 

BC  Downstream  control  edge 

BSP  Splitter  width 

c  Control 

cl  Center  line 

core  Pertaining  to  potential  core 

d  Delay 

e  Entrainment 

eff  Effective 

enet  Net  entrainment 

j  Jet 

L  Loop 

n  Nozzle 

o  Oitput  pressure  and  width  for  BQ,  PQ;  orifice  for  all  others 

r  Relative 
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NOMENCLATURE  (Cont'd) 


rec  Recovered  (block  loaded) 

s  Supnly 

sp  Splitter 

SB  Spill  back 

t  Throat 

th  Nozzle  throat 

T  Total 

v  Vent 

vo  Virtual  origin 

XSP  At  splitter 

y  Deflection 

• 

0  Due  to  rotation 

0  Reference  0 

1  Reference  1 

2  Reference  2 

3  Reference  3 

90  -90  deg  of  phase  lag 

Superscripts 

—  Average 

’  Dimensional  quantities  in  temperature  compensation  equations 
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APPENDIX  A.— COMPUTER  PROGRAMS  FOR  LJARS 


This  appendix  lists  computer  programs  for  the  laminar  jet  angular 
rate  sensor  (LJARS)  (listings  A-1  to  A-3),  defines  the  terms  (table 
A-1),  and  shows  a  sample  printout  (fig.  A-1)  and  a  standard  silhouette 
(fig.  A-2)  of  the  laminar  proportional  amplifier. 


APPENDIX  A 


LISTING  A-1.  LJARS 


10PRINT  "IF  YOU  ARE  INTERESTED  ONLY  IN  AN  ANGULAR  RATE  SENSOR  AN 
SVER  NO  FOR  PLOTS l INPUT  "TODAYS  DATE  , C$ 

20SELECT  PLOT  413. PRINT  213 (100 )t PLOT  [,,]:SELECT  CO  213(100) 
30PRINT  C$ (PRINT  "LAMINAR  JET  DEVICE  CHARACTERISTICS" 

40DIM  V5(56),W5(56),V6(56),W6(56),V7(56),W7(56),C7(55),G1(100),P 
4(100), C(25),FO(25) 

50 INPUT  ADO  YOU  WANT  ANGULAR  RATE  DATA" ,R$ 

60 INPUT  "DOES  THIS  DEVICE  HAVE  CENTER  DUMP  OR  OUTPUT  VENTS", Y$:I 
P  Y$“"YES"THEN  70 : V1-1E-10 : V2- 1E10 :L8-1E10 :L 9- 0 

70Y2**0  SlNPUT  "DOES  THIS  DEVICE  HAVE  SUDDEN  EXPANSION  OUTPUTS",  Y1 
$:IF  Y1$M"YES"THEN  80sY2-l 

80 INPUT  "MODIPIED(l)  OR  WIDTH(2)  REYNOLDS  NO.",Y3 
90 INPUT  "DO  YOU  WANT  RANGES, DO  YOU  WANT  PLOTS" , Z 3$, Z 9$ 

100 IF  Z3$-"YES"THEN  270:IF  Z9$-"NO"THEN  120 

1101NPUT  "RESPONSE  OTHER  THAN  SELF-STAGED" , Z 5$ : IF  Z5$-"NO"THEN  1 
20:INPUT  "BIAS  PRESSURE  FOR  RE SPONSE" , J 6 s INPUT  "SELF-LOAD ( 1 ), BLO 
CKED (2 ) , BOTH ( 3 )" ,Z2 

120 INPUT  "CONTROL  AND  THROAT  WIDTH", W,B 
130IF  Y3-1THEN  160 

140 INPUT  "ASPECT  RATIO  AND  REY.NO.",S,N 
150GOTO  170 

160INPUT  "ASPECT  RATIO, MOD  REY . NO . " , S ,N9 
170INPUT  "SPLITTER  DISTANCE  ",B6 

180INPUT  "NOZZLE  AND  SPLITTER  LGTH",L4,C9 
190 INPUT  "OUTPUT  LENGTH , AVG . WIDTH  ",X3,F8 
200IF  Y$ l 1 "YES"THEN  220 

210INPUT  *CNTR  DUMP  W.L;  VENT  W, L  ", VI, V2, L8, L9 
22 0 INPUT  "BIAS  CONTROL  PRESSURE  ",H3 
230INPUT  "CTRL  L , AVG . W, MIN . W, ENT . R" , XI , Wo, W3, B1 
240 INPUT  "SPLITTER  ,  OUTPUT  WIDTH  " , B5, E7 :PRINT  "  " 

250IF  Z3$-"YES"THEN  310 
260G0T0  390 

270INPUT  "RANGE  AND  STEP  OF  NR',BC,BT,AR,XSP",N1,N2,N3,Z4,Z5,Z6, 

Z7,Z8,Z9,S3,S4,S5,X4,X5,X6 

280J-0 

290GOTO  180 

300REM  LOOPS  FOR  RANGES 
310FOR  N9-N1T0  N2STEP  N3 
320FOR  S-S3TO  S4 STEP  S5 
330F0R  B6-X4TO  X5STEP  X6 
340F0R  W-Z4TO  Z5STZP  Z6 
350FOR  B-Z7TO  Z 8 STEP  Z9 
36 OIF  J-OTHEN  390 

370SELECT  PLOT  413:PLOT  ( , -3 95  , ] : SELECT  PRINT  213 
380REM  CALL  DISCHARGE  COEFF.  ROUTINE 

390J-1:IF  Y3-1THEN  400 : N9-N / (L 4+1 ) / (1+1/ S) I  2 : GOSUB  4020:GOTO  410 


400G0SUB  4020:N-N9*(L4+1)*(1+1/S) 12 

410IF  N*S  [ 1400THEN  430:PRINT  "****  WARNING  ****":PRINT  "MA 
X.DYN. RANGE  REYNOLDS  NO . EXCEEDED  ...  FLOW  BECOMING  TURBULENT" 
420REM  COMPUTE  2-D  CD , ENTR. FLOW, JET-EDGE  SPACE , SP ILLBACK  FLOW,  N 
ET  ENTR. COEFF ., SPILLBACK  LOSS  COEFF. 

430D“ (1- S) /2+SQR (( S  +  l ) I 2/4-S*(l-C)  ) 

440Q-1 ,651*(.02772+1.32*(W+B1)/(C*N))I (1/3)-. 5 

450K-1/ (C*Q+D/2) :Y-K*B/2:B8».5*(B-(1+W/.021/C/N) 1(2/3)) 

460T1"LOG((1+.5*D*K)/ ( 1- . 5 *D*K) )  /  2 

470K3-4/ (EXP (T 1 )+EXP (-T1 ) ) I  2 : K4-K 3/ Q/C : R4- 0 

48  0Q4- Q*EXP(-K4*B8)/2 ! A-K4 *Q4* 2 : C 3- 1-3 . 0303+K4 *Q4 1  2 

490REM  IF  CONTROL  PARAMETERS  UNCHANGED  BYPASS  RV  CALC. 

500IF  W[ 1 M5THEN  560 
510IF  B  1M6THEN  560 
52 OIF  N  1 M7THEN  560 
530 IF  S  JM8THEN  560 
540 GOTO  660 

550REM  COMPUTE  VENT  RESISTANCE 
560FOR  P1-.1TO  , 5STEP  .05 

570P-P11 .25 ! S 1-LOG (1/P+SQR(1/P! 2-1)) S Y 0- S 1/ Ks D 1- (B /2 -Y0) /I 0 s 1-0 
580FOR  Y1-Y0T0  B/2 STEP  Dl 
590S2-2 /  (EXP(Y1*K)+EXP(-Y1*K)) 

600A1- (SQR(S2! 4+Pl)-S2! 2)*D1 :I-I+A1:NEXT  Y1 
610R1-D+P1/I 

6202RV-  t.ttt,  NR-  tttt.t,  NR  tt  ti.t  ,AR-  « .  «*l  9  ,  BC-  ##  .  ,  BT-  Iff  .  ##  , 

XSP-#d.0C,CD-  .tUt 

630R4- R4+R1 

640NEXT  PI 

650R1-R4/10 
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LISTING  A-1.  LJARS  (Cont'd) 

660M5-WsM6-BsM7-N:M8-S 

67 OREM  CHANNEL  RESISTANCE  CORRECTION  AND  COEFFICIENTS 

68 OIF  W0/S]2THEN  690 iK2- . 35 sGOTO  700 

690K2-.5 

700T4-24*X1*C*(S/W0+W0/S+K2)/(W012*S*N) sT5-. 95*(C/W3) 12  SC2-1.32* 

c*c 

710IF  V1/S12THEN  750 
720IF  S/Vl] 2THEN  750 
730K9-. 35 
740GOTO  830 
750K9-.5 

760  IF  C ] . 5THEN  770 : N0-2*C : GOTO  790 
770  NO-1 

780REM  JET  CHARACTERISTICS , D I SP . TU , E FF  .  NOZZLE  L. , CORE  L., VIRTUAL 
ORIGIN, MOMENTUM  LOSS  COEFF . , RATE  SENSOR  EFF . L . , PART ICLE  TRANSIT 
TIME 

790  DO- ( 1+S) /4-.5*SQR((l+S)12/4-(l-C)*S) 

800  X-N*(D0*10/3) 12/25 *N0 
810  XO-X*  ( •  2  5  / (DO*1.4)!2-l) 

820  X2-. 0155*C*N SPRINT  "CORE  LENGTH-" , XO, "V IRTUAL  ORIGIN-" , X2 , "D 
ISPLACEMENT  THIC KNE S S-" , DO 

825  IF  S [ 1THEN  82 6 : N 5- 4 : N6- 5 . 6 sGOTO  830 

826  N5-3.00:N6-3.00 
830IF  XO] « OTHEN  850 

840  CO- (1+B6/X2) ! (-2 / 3 ) s NO- 2 *C s Z 0- 9*X2 I  2 * ( ( 1+B6/ X2 ) !  (7/3)-l)/28- 
. 75 *B6*X2 s FO- 1 /I 2 /X2 / ((1+B6/X2)!  (4/ 3 ) -1 ) s GOTO  90  0 

850  IF  X0]B6  THEN  880sIF  C ] . 5THEN  86 0 s C 0- ( 1+ ( B6-X0) /X2 ) I (-2 /3 ) s G 
OTO  870 

860  CO- (1+(B6-X0)/X2) ! (-2/ 3 ) * ( 1-N5*SQR ( (X+XO) /N) / S) / ( 1-N6*SQR (X/ 
N)  /S) 

87  0  Z0-X0*(B6-X0/2)-. 7 5 *X2 * <B6-X0)+ 9*X2 ! 2 * ( ( l  +  ( B6-X0) /X2 ) 1(7/3)- 

1) /28sF0-l/16/(X0+.75*X2*((l+(B6-X0)/X2)I  (4/ 3 ) -1 ) ) 5 GOTO  900 
880C0- (1-N5*SQR( (B  6+X )/N)/S)/ (1-N6*SQR(X/N)/S) sNO-l s ZO- B6 1  2 / 2 s FO 
-1/16/B6 

890REM  PRESSURE  RECOVERY  AT  +-1Z  AND  0 

900P2-K3sG0SUB  4330sG0SUB  4280 s A3-Q2+. 5 

910B2-2/  (1+B6/.021/C/N)!  (2/3)sIF  B7[B2  THEN  920sB2«B7 

920P6-C0*(1-1.1*B5/(SQR(C*N*B5/2)))*C2/B2*(1-SQR(C2*C0)*V1I2*S*N 

/(24*V2*(S/V1+V1/S+K9)))*(1+2*C*P1*SQR(ABS(P1))/  (R1*C2))*C3 

930D6“C0*(1-1.1*B5/(SQR(C*N*B5/2)))*W*W/B2*(1-SQR(C2*C0)*V1I2*S* 

N/(24*V2*(S/V1+V1/S+K9)))*C3 

940D5«24*L9*C*(S/L8+L8/S+.5)/S/N/L8! 2 s J 9-1+ . 5  * (D5+. 95*C*C* . 5 /L8 ! 

2) /P6 

95 OGOSUB  4330  sGOSUB  4190 
960E-1 S02-1 S03-1 
9 7 0 RE M  OUTPUT  CHARACTERISTICS 
980F0R  M-~. 01  TO  .01STEP  .01 
990 IF  Z9$-"YES"THEN  1000sM-0 
1000P5-P6+G*M/2 
1010A5-A3*P5/P6 

1020R5»24*X3*C*(S/F8+F8/S+K2)/(F8! 2*S*N) 

1030K5-. 95  *C ! 2 / B  7 ! 2 
1040IF  M [ ] 0 THEN  1090 

1050F0R  I-1T0  5  s  PRINT  HEX(OA)sNEXT  I 

1060PRINT  "  CENTERED  INPUT  AND  OUTPUT  CHARACTERISTICS" 

1070  FOR  I-1T0  3  s  PRINT  HEX(OA)sNEXT  I 

1080PRINT  "  Q  PO  PO *Q0  RO  PC  RC" 

1090H-0 

1100E2-0 sE3-0sE4-0sE6-0sE7-0:E8-0:T7-0sA6-0sZ3-l 
11 10 FOR  QO-OTO  . 5STEP  .OlOsIF  A6-1THEN  1170 
1120R-D5+.95*C*C*( .5-QO) /L8I 2 
1130IF  S [ 2THEN  1140:Z3-2/S 

114OPO-P5*(l-2*QO*C/SQR(P5*B7))/(l-Q0*C/B7/SQR(B7*P5))-R5*QO-K5* 
Q012+C*SQR(P5/B7)*Q0*.25*S*Z3+(.5-Q0)*R+(C*00/B7) 12*4*Y2 
1150R6-C*SQR(P5/B7) / ( 1-Q0*C / B7 / SQR (P5 *B 7 ) ) * ( -2+ ( 1-2 *QO*C / SQR (P 5* 
B7)) / B7/ (1-Q0*C/B7/SQR(P5*B7)) ) -R5- 2*K5*Q0+C* SQR (P 5/ B7 )* . 2 5 *S*Z3 
-R-. 95*C*C*(.5-Q0) /L81 2+8*(C/B7) I 2*Q0*Y2 

1160Z-R5+K5*Q  0+P  0/ (Q0+. 000001) sA4“ . 5/A5sIF  A4 [ZTHEN  1180sIF  P0/( 
A5-Q0) [ -R6THEN  1180  s A6-1 s A7-Q0 5 A8-P 0 
1170P0-A8/ (A5-A7)*(A5-Q0)sR6--A8/ (A5-A7) 

1180E1- (P0-E2)*1000+E3 
1190E2-P0 

1200E3-E1-INT(ABS(E1))*SGN(E1) 

1210E4-E4+INT(ABS(E1))*SGN(E1) 

1220E5- (Q0-E6) *500+E  7 
1230E6-Q0 
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LISTING  A-1.  LJARS  (Cont'd) 


f 

\ 

k- 


i 

[.  i 


1240E7“E5-1NT(ABS(E5))*SGN(E5) 

1250E8-E8+INT(ABS(E5))*SGN(E5) 

17  60REM'  STORE  DATA  POINTS  FOR  OUTPUT  PLOTS 


12 70 IP  M]-.'01THEN  12  80  s  V5  (E  )-El  s  W5  (E)-E5  s  L5-E4  :  Y5-E8  :  E-E+l :  GOTO 
1300 

1280IF  M]0THEN  1290 :V6 (02 )-El :W6 (02 )-E5 :L6-E4 :Y 6-E8 : 02-02+1 : GOTO 
1300 

1290V7 (03)-El:W7(03)-E5:L7“E4:Y7-E8  s 03-03+1 

1300IP  M[ J  OTHEN  1450 

1310Q2-Q0 

1320G0SUB  4000 

1330G0SUB  4260 

1340P-P0-P2 

1350IP  P ] OTHEN  1420 

1360IP  T 7-1 THEN  142 0 :T7-1 sR7- R6 

1370P7-P0 

1380R8-R0 

1390Q7-Q2 

1400  PRINTUSING  1410, QO 

1410Z  SELF-STAGED  OPERATING  FLOW  -  -t.ttt 
1420IF  (H-INT(H) ) ]0THEN  1450 
1430PRINTUSING  144 0 , Q 0 , P 0, QO*P 0, R6 , P2 , RO 
1440 Z-t.tt  -#.####  -#.##«#-«.###  -#.### 

1450H-H+.2 

1460IF  PO [ OTHEN  1480 
147 ONE XT  QO 

148 OIF  Z  9$-"YES" THEN  1490:M-.01 
1490NEXT  M 

1500IF  Z9$-"N0"THEN  1760 
1510REM  DRAW  AXES  AND  PLOT  OUTPUT 
1520SELECT  PLOT  413 

1530PLOT  [500,300, 1,60[,-5,HEX(27) ],100[5,,HEX(2D) ], [,,"+"], 5 [ -1 
00, ,  "+"  ] ,  3  [ ,  100,  "+*’  ]  ,  [  ,-300,"+"  ] 

1540FOR  E9-1T0  E-l 
1550PL0T  [V5(E9),W5(E9),"*"] 

15 6 ONE XT  E  9 

1570PLOT  [ -L5 , -Y  5, "+" ] 

1580FOR.  E9-1T0  02-l:PL0T  [ V 6 (E9) ,W 6 (E9) , ] :NEXT  E9 
1590PLOT  [ -L6 , -Y  6 , "+" ] 

1600F0R  E 9-1 TO  03-l:PL0T  (V7 (E9) , W7 (E9) , ] :NEXT  E9 
1610PLOT  [ -L7 , -Y7 , "+" ] 

1620REM  PLOT  INPUT 
163004-0:02-0:03-0 


1640C4-INT(100*Q7) /100+.04 
1650F0R  Q2-OTO  C4STEP  .002 
1660GOSUB  4000 
1670 GO SUB  4260 
168001- (P2-02)*1000+03 
169  002” P2 

170003-0 1-INT(ABS (01 ))*SGH (01) 

171004-04+INT (ABS (01 ) ) *SGN (01) 

1720PLOT  [01,1,"+''] 

17 30NEXT  Q2 

1740PLOT  [ -04, -(C 4/.  002+1 ),"+"] 

1750REM  COMPUTE  DEFL. RESIST. 

1760R0-R8 
1770P2-P7 
1780Q2-Q7 
1790P1-P2-R0*Q2 
1800  GOSUB  4190 
1810  GOSUB  4130 
1820G2-G*R3/ (R3-R7) 

1830IF  Z  9$-"N0"THEN  1930 
1840REM  PLOT  DEFL. RESIST. 

1850PLOT  [P7*1000-R3*150,Q7*500-75,HEX(FB)1:I4“0 
1860FOR  I1-1T0  30:PLOT  [ R3*10+I 3 , 5 , HEX (FB ) ] : I 2-R2 
T(I2 ) : 14- I4+INT( 12 ) 

1870NEXT  I1:PL0T  [ -P 7*1000+R3*15 0-14 ,-Q7*500-75 , "+" ] 

1880REM  LABEL  AXES 

1890PLOT  [ 0, -15 , "0" 1 , [ 80 , 0 , "0. 1"  1 , [70,0,"0.2"1 
70,0,ft0. 3"], [70,0, "0.4"], [70,0," 0.5"] 

-545,15," 0"), [-10, 100, "0.2"] 
i-30, 100,  "0.4"], [-30, 100,  "0.6"],  [  205  , -330 ,  "PRESSURE"  ]  ,  [ 
-335, 305, ftFL0W"] , [-250,-325, ] 

1930SELECT  PRINT  213:PRINT  HEX ( OA) : PR INTUS ING  62 0 , R1 , N, N9 , S , W, B , 
B6,C 


13*10+13: I3-I2-IN 


1900PLOT 

1910PL0T 

1920PL0T 


LISTING  A-1.  LJARS  (Cont'd) 
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1940PRINT  *'  " 

1950REM  EFFECT  OF  CONTROLS  ON, AND  RATE  SENSITIVITIES 
I960  ZO-ZO+W*W/2*(2/C*SQR(RO*Rl/ (R0+R1) ) * (Q-Q4-H3/R0) /SQR(ABS (Q- 
Q4-H3/R0)) -A) *(2*B6/W-l)/(4  *C2*( R0+R1) /W/W/R0/Rl-(2/C*SQR(R0*Rl/ 
Q4-H3/R0)/ SQR (ABS (Q-Q4-H3/R0)) -A) ) 
4*P6*Z0/N0/N/C:S8-2*N*P6*Z0/C/57. 3/NO 

2090,"  SELF-STAGED  DEFL.  RESIST. -",R3 
2090, "SELF-STAGED  OUTPUT  RESIST. -",R7 
2090, "BLOCKED  GAIN  AT  SELF-BIAS  -”,G 
2090,"  SELF-STAGED  BIAS  -”,P2 

2090,"  SELF-STAGED  GAIN  -H,G2 

2030IF  R$-"NO"THEN  2110 
2040PRINTUSING  2100, "NORMLIZED  RATE 

- -  RATE 

RATE 
RATE 


(R0+R1) ) * (Q* 
1970  S 7* 
1980PRINTUSING 
199 OPRINTU SING 
2000 PRINTUSING 
201 OPRINTU SING 
202 OPRINTU SING 


2100, "DIMNSLESS 
2100, "MILH  5606 
2100, "PNEUMATIC 


SENSITIVITY-", S7 
SENSITIVITY-", S8 
SENSITIVITY-", S8* 2. 283  5 E-6 
SENSITIVITY-", S8*2.649E-9 


2050PRINTUSING 
20 60 PRINTUSING 

207  OPRINTU SING _ 

2 08 OPRINTU SING  2100, "SENSR-90DEG  FREQ , F*BS I 2/NU-" , FO*C*N : IF  Z9$- 
"NO"THEN  2960 

2090*###########  tit  tiiit  #########  -###.### 

2100 *############################  -#.### Ill  I 
2110PRINT 


BIAS  PRESS. 


P  JET 


GAIN 


DEFL. RES." 


2120PRINT 
2130G6-0 

2140REM  COMPUTE  GAIN  VS. BIAS  AND  STORE  DATA 

2150G4-0 :G8-1 :FOR  P2— .05T0  . 1 95STEP  . 005 :GOSUB  4330:GOSUB  4280: 
COSUB  4190 :G7  (G8  )-G:  G8-G8+1:  IF  ABS  (G  )  [  99TI1EN  2160:G6-100 
2160NEXT  P2 

2170F0R  G8-1T0  49STEP  10: P2- (G8-11) * . 005 : GOSUB  4330:GOSUB  4280:G 
-G7 (G8) : GOSUB  4130:COSUB  4300:NEXT  G8 

2180IF  G 6 [ 100THEN  2190:S£LECT  PLOT  413:PLOT  [ 80, -125 ,]: SELECT  PR 
INT  213: GOTO  2280 
2190IF  Z9$-"N0"THEN  2960 
2200REM  DRAW  AXES  FOR  G  VS.PCBIAS 
2210SELECT  PLOT  413: 

2220PLOT 
22  30PL0T 

2240REM  PLOT  G  VS.PCBIAS 
2250G9-0:G4-0:FOR  G8-1T0  50:J7- (G7 (G8)-G4)*10+G5 :G5-J7-INT(ABS(J 
7))*SGN(J7):G 9- G 9+1 NT (ABS(J7))*SGN(J7):PL0T  [ 10 , J7 ,"*"]: G4-G7 (G8 


_  _ :PLOT  [120,0,] ,  50[ 0, -5 ,HEX (27 ) 1 

80[5,0,HEX(2D)] , [ , ,"+"].5[-100,0,"+"] ,20[5,0,HEX(2D) ] 

[0,2  00,"+"],2[0,-i00,"+,,].  [-110,,  ] 


20, -15, "-.05"], [80,0," 0"], [80 
, [70,0, ".20"] 

20" ], [70, -230, "CONTROL  BIAS  P 


2300 IF  Z5$-"N0"THEN  2310 : P2-J6 :G0 SUB  43 30 : Q7-Q2 : P 7- J6 :  IF  Z2-1THE 
N  2310  :  IF  Z2-3THEN  2 310 : X3- 0 : R7-0 
2310C5-WI3/ (12*C2) 

2320L1-2*C*C* (X1/W0+X3/F8) 

2330Q2-Q7 
2340G0SUB  4000 
2350R9-R0-R7 

2360K6-W*W/(4*C2):P1-P7-R0*Q2 

237  OKI- (1+R9*(A*K6+1/R1+2*K6*P1/ ( S QR (ABS (PI ) ) *C ) ) ) / (2*L1*C5 ) 

2380 AO- K1*2*L 1*C  5 
23 90 A3- (A0-1)*L1/R9 
2400A2-2*L1*C5 

241OKO-1+R0* ( A*K6+1/ R1+2*K6*P1/ (SQR (ABS (PI) )*C) ) 

2420G3-1 :P  3-0 
2430Z3-0 
2440G0SUB  4190 
2450T8-0 

2460FOR  L-1T0  100 
24  70F2-10 I  (L/50-3) 

2480W9-2*#PI*F2 

2490G0-G*K0/SQR((A0-A2*W9! 2)12+(A3*W9) 12) 

2500F1-SQR(K1) / (2*#PI) 

2510G1 (L ) -20*. 4 32 94 5* (LOG (GO) -LOG (G3)) 

252  0C3-G0 

2530P8-ARCTAN(A3/SQR(A2*A0)*F2/F1/ (1— (F2/F1)  I  2)) 

2540IF  P 8] OTHEN  2560 
2550Z3-180 

2560P9-P8*180/#PI+B6*F2*360+Z3 

2570P4(L)-P9-P3:P3-P9 

2580IF  P  9-45 [ OTHEN  2600 


2260NEXT  G8:PL0T  [-3 90 , -G9, "+" ] , [ -1 
,0.".05"], [70,0, ".10"], [70,0, ".15"! 
2270PL0T  [-4 50, 115, "10"], [-20,100, 
RESSIIRE"],  [-340, 200,  "GAIN'1] 

[J"v -  ' 


RESSURE 
2280 IF  Z  3$ 


'YES"THEN  2290: Z5$-"N0" 


2290SELECT  PLOT  41 3 : T6-0 : Z6S-" 1" 
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LISTING  A-1.  LJARS  (Cont'd) 


2590IF  T8-1THEN  2600:T8-1:F3-F2 SG2-G0 
2600NEXT  L 

2610IF  T6-0THEN  2620 ! Z1-G2+F3 s GOTO  2670 

2620Z«G2*F3 

2630SELECT  PLOT  413 

2640PLOT  [500,230,] ,80 [0,-5, HEX (27)], 100t 5, 0, HEX (2D) ] 

2650PLOT  80 [ 0 , 5 , HEX(27 ) ] 

2660PL0T  [0, 0, "+"], 4[0, -100,"+ "],2[-250,0,"+"], [120, ,"+"], [250,, 
"+"], C-370,,],4[0,100,"+"] 

2670PLOT  [0,-400,+']:U2“0:U3"0 
2680F0R  L-1T0  100 
2690U1-G1(L)*10+U2 
2700U2-U1-INT(ABS(U1))*SGN(U1) 

2710U3-U3+INT(ABS(U1))*SGN(U1) 

2720PLOT  [ 5 , U1 , Z  6$ ] 

2730NEXT  L 

2740PLOT  [-500,-U3,"+"]:U5»0:U6“0:PLOT  [0,400,"+"] 

2750FOR  L-1T0  100 
2760U4«P4(L)*2+U5 
2770U5»U4-INT(ABS(U4))*SGN(U4) 

2780IF  U 6  [400THEN  2800 
2790G0TO  2830 

2800U6-U6+INT(ABS(U4)) *SGN (04) 

2810PL0T  [5,-U4,Z6$] 

2820NEXT  L 

2830PL0T  [-5ML-1)  ,06,"+"]  :  IF  Z2-3THEN  2840sPL0T  [ -40 , 0 , "40" ] , [ - 
20, -100, "30"] :GOTO  2850 
2840Z2-2sT6-l:Z6$-"2"  -.GOTO  2300 

2850PLOT  [-20, -100, "20"] , [-20.-100, "10"] , [ -10, -100, "0" ] 

2860PL0T  [0,-20,".001"],[85,,.003" ],[95,0," .01"],[85,,".03"],[1 
05 , 0 , "  .  1 "  J 


28  70PL0T 
2880PLOT 
2890PLOT 


"-200"1 , [-40,100, "-150"] , M 
-40 , 100, "-50" ] , [ -20, 1 00 , "0" ] , [ ,-150; "PHASE" ] 


10,20,"-: 


[-40, 100," -100"] 


-440,-290, "NORMALIZED  FREQUENCY"] 

2 9 00 IF  Z 5$”" YES" THEN  2910:PLOT  [ -4 50 , 2 50 , "SELF" ], [-40, -15 , "STAGE 
D"j, [-60,-15, "GAIN  DB" ] : GOTO  2920 

2910PLOT  [-450, 250, "DB  GAIN" ], [-30 , -3 0 ,],[ 100, -145 ,] s SELECT  PRIN 
T  213 : GOTO  2930 

2920PL0T  [100,-145,] jSELECT  PRINT  213 
2930IF  Z-OTHEN  2 94 0 : T 7-1 : PRI NTU S INC  2950, T7.Z 

2940IF  Z1-0THEN  2960 s T7- 2 : S ELECT  PLOT  413:PLOT  [ 600 ,,]: SELECT  PR 
INT  213:PRINTUSING  2950,T7,Z1 
2950ZGBP  (#)  —#.##### 

2960IF  Z 3$“ "YES" THEN  3020 

297 0PR1NT  "  ":PRINT  ”  ":PRINT  "  ":PRINT  "  "sPRINT  "  " 

2980INPUT  "DO  YOU  WANT  TO  DO  ANOTHER  PROBLEM", Z4S 
2990IF  Z4$-"NO"THEN  3040 

3000SELECT  PLOT  413:PLOT  [ ,-95 ,]: SELECT  PRINT  213 
3010GOTO  30 

3020NEXT  B : NEXT  W:NEaT  B6:NEXT  S:NEXT  N9 
3030PRINT  "  " : PRINT  "  "sPRINT  "  " 

3040 SELECT  CO  005 (64 ), PRINT  005(64):END 

4000R0-T4+T5*Q2 

4010RETURN 

4020  FOR  I0-1T0  25 

4030  IF  I0]1  THEN  4050 

4040  C(I0)-.5 :G0T0  4070 

4050  IF  1 0] 2  THEN  4070 

4060  C(I0)-.6 

4070  F0(I0)-2.667/C(I0)*(EXP (6. 6*C (10) I 2)-l)-N9 
4080  IF  I  0 [ 2  THEN  4110 

4090  C (1 0+1 )“C(I0)-F0(I0)*(C(I0)-C(I0-1))/(F0(I0)-F0(I0-1)) 

4100  IF  ABS (FO (1 0) ) [ , 01  THEN  4120 
4110  NEXT  10 
4120  C “C ( 1 0) : RETURN 
4130  K7«(W*W)/(4*C2) 

4140  IF  Q2J0THEN  4 15 0 : B2-0 : GOTO  4160 

4150  B2-3.4242*(l/ (1+B/ 0. 1 /Q2/C/N ) -0. 5 )  :  IF  B 2 ] OTHEN  4160:B2-0 
4160  D2"(l-((B2*Q2*A)/2)-((B2*Q2*Pl)/ (SQR( ABS (PI ) )*C ) ) ) / (2*K7+(B 
2*Q2)/R1) 

4170R3-(D2*(1+R0/R1)+(R0*P1)/(C*SQR(ABS (Pl)))+R0*A/2)/(D2*(l/Rl) 
+A/2+(Pl/ (SQR(ABS(P1))*C))) 

4180RETURN 

4190  K7-(W*W)/  (4*C2) 

4200  IF  Q2J0THEN  4 2 1 0 : B2-0 : GOTO  4220 

4210  B2-3.4242*(1/(1+B/0.1/Q2/C/N)-0.5):IF  B2] OTHEN  4220tB2-0 


110 
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LISTING  A-1.  LJARS  (Cont'd) 

4220  D2«(l-((B2*Q2*A)/2)-((B2*Q2*Pl)/(SQR(ABS(Pl))*C)))/(2*K7+(B 
2*Q2)/R1) 

4230  J3-IsIF  Q21-0THEN  4240sJ3-0 

4240G-D6* (1/ (D2*2* ( 1+R0/ R1 )+R0*A+(R0*  2*P1 ) / (C* SQR (ABS (PI) ) ) ) ) * J 9 
*1/K7*(2*B6/W/(1+2*Q2*J3)-1) 

4250RETURN 

4260P2—  (Q-Q4)*R1+(R1+R0)*Q2 
4270RETURN 

4280P1-  (P2/R0-(Q-Q4) )*R0*R1/ (R0+R1) 

4290RETURN 

4300PRINTUSING  431 0 ,P 2 , PI , G, R3 
4310Z  -#.###  -i.titli  -###.## 

4320RETURN 


4J20UTURII 

4330Q2— .5*(T4+R1)/T5+SQR(.25*((T4+R1)/T5) ! 2+ (P2+R1* ( Q-Q4 ) ) /T5) : 

GOSUB  4000 

4340RETURN 


LISTING  A-2.  COMPUTER  PROGRAM  FOR  SIMPLIFIED  LuARS  SENSITIVITY 

5  SELECT  PLOT  413, PRINT  213 

6  PLOT  (,,  ] 

7  SELECT  CO  213 

10  DIM  C ( 25 ) , F ( 25 ) 

20  INPUT  "ASPECT  RATIO", A 
30  INPUT  "NOZZLE  LENGTH" ,L 

35  INPUT  "NOZZLE  WIDTH  IN  MM" , B3: B 3- B3/ 1 000 
40  INPUT  "OUTPUT  WIDTH", BO 
50  INPUT  "SPLITTER  WIDTH", B1 
60  INPUT  "SPLITTER  DISTANCE" ,X 

70  INPUT  "INPUT  WORKING  FLUID" ,A $ J IF  A$-"AIR"  THEN  80:IF  A$«"OIL 
"  THEN  90 

80  M-1.789E-5 :R-1 . 2 059 :GOTO  100 
90  M-. 0155 :R“860 

100  PRINT  "  MOD. NR  NR  SIGMA*NR  S(PA/0/S)  S(PSI/0/S)  P 

RES  SURE (PA)" 

110  FOR  N-50  TO  150  STEP  10 
120  Nl-N*(l+L)*(l+1/A)  !2 
130  GOSUB  1000 
140  C2“l . 32  *C 1 *C 1 

150  S-M*N1*C2*X*X* (1-8*X/C2/N1/A/A)*(1-1. 1 *B 1/ SOR (Cl *N 1*B1/ 2 ) ) /5 
7.3/ Cl/ BO 

160  Sl-S/1000/6.89 

165  P«R/2*(N1*M/R/B3) !2 

170  PRINTUSING  18 0 , N , N1 , A*H1 , S , S 1,  P 

180  2  HOtt.H  ttttf.it  Iff.  if  9*9  #.#<•#!!!!  ttttftt . 

ttt 

190  NEXT  N 

195  SELECT  CO  005, PRINT  005 
200  F.ND 

1000  FOR  I-1TO  25 
1010  IF  IJ1THEN  1030 
1020  C (I ) “. 5 (GOTO  1050 
1030  IF  I ] 2 THEN  1050 
1040  C( I ) -. 6 

1050  F(I)-2.667/C(I)*(EXP(6.6*C(I) I2)-l)-N 
1060  IF  I [ 2THEN  1090 

1070  C(I+1)-C(I)-F(I)* (C(I)-C(I-l )) / (F (I)-F(I-l) ) 

1080  IF  ABS(F(I)) [.01THEN  1100 
1090  NEXT  I 

11  00  Cl-C(I) : RETURN 
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LISTING  A-2.  COMPUTER  PROGRAM  FOR  SIMPLIFIED  LJARS  SENSITIVITY  (Cont'd) 

ASPECT  RATIO?  1.6667 
NOZZLE  LENGTH?  2 
NOZZLE  WIDTH  IN  MM?  .75 
OUTPUT  WIDTH?  1.7 
SPLITTER  WIDTH?  .5 
SPLITTER  DISTANCE?  20 
INPUT  WORKING  FLUID?  AIR 


MOD. NR 

NR 

SIGMA*NR 

S(PA/G/S) 

S(PSI/0/S) 

PRESSURE (PA) 

50 

383. 99 

640. 00 

0.  01509 

2. 1 9  OE -06 

34.7  8  6 

60 

460.79 

768. 00 

0.  02  03  0 

3.  020E-06 

50.092 

70 

537. 59 

896. 00 

0.  02662 

3.  864E-06 

68.180 

30 

614 . 3  9 

1024.00 

0.  03254 

4.  723E-06 

89. 052 

90 

691.18 

1152. 00 

0.  03  3  54 

5.  594E-06 

112.707 

100 

767. 98 

1280. 00 

0.  04462 

6.476E-06 

139.144 

no 

844  .  78 

1403. 00 

0.  05  07  6 

7.  363E-06 

168. 364 

120 

921. 58 

1536.00 

0.  05  C  97 

8.  269E-06 

200.  368 

130 

998. 38 

1664.00 

0. 06324 

9. 178E-9C 

235. 154 

14  0 

1075.18 

1792.00 

0.  06956 

1. 009E-05 

272.723 

150 

1151. 98 

1920. 00 

0. 07592 

1. 1G2E-05 

313. 075 

ASPECT  RATIO?  2.5 
NOZZLE  LENGTH?  2 
NOZZLE  WIDTH  IN  MM?  .75 
OUTPUT  WIDTH?  1.7 
SPLITTER  WIDTH?  .25 


SPLITTER  DISTANCE?  20 
INPUT  WORKING  FLUID?  OIL 
vnn  ”B  "P  crcMi* 


MOD. NR 

NR 

SIGMA*NR 

S(PA/0/S) 

SCPSI/O/S) 

PRESSURE (PA) 

50 

294.00 

735. 00 

11.94935 

1.  734E-C3 

21463.823 

60 

352.80 

882. 00 

15.71710 

2.  281E-03 

30907. 905 

70 

411. 60 

1029.00 

19.56259 

2. 839E-03 

42069.093 

30 

470.40 

1176.00 

23.47678 

3.407E-03 

54947.387 

90 

529.20 

1323.00 

27.45065 

3.  984  E -03 

69542. 787 

100 

588. 00 

1470.00 

31.47818 

4.  568 E -03 

85855. 293 

110 

646. 80 

1617.00 

35.  55275 

5.  160E-03 

103884. 904 

12  0 

705. 60 

1764.00 

39.67135 

5.  757E-03 

123631. 621 

130 

764. 40 

1911.00 

43.82959 

6.  361E-03 

145095.445 

14  0 

82  3.2  0 

2058.00 

48.02430 

6.  970E-03 

163276.374 

150 

882.00 

2205. 00 

52. 25240 

7.  583E-03 

193174.409 
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LISTING  A-3.  COMPUTER  PROGRAM  FOR  REYNOLDS  NUMBER  CONTROL 


5  REM  *******  SHUNT  RESISTANCE  COliPLUSAT IOI1  ******* 

10  DIM  C(25),F(25) 

30  INPUT  "NOMINAL  MODIFIED  REYNOLDS  NO."  N 
40  INPUT  "VISCOSITY  CHANGE  RATIO", V 

SO  PRINT  "NOMINAL  R (EOZ ) /R(C AP )  NR*  AT  NEW  VISCOSITY  NR/NRO 

60  FOR  R'.  0001  TO  5.0001  STET  .50 
70  N9«N:COSUB  4010 
80  C0-C1:C-C1 

90^N9-2.667*(EXP(6.6*C0*CO)-l)  /C 0/CO*C  •'■  ( SQR  ( .  25+V*R*  (R+l ) *C0 1  2 /C 

100  GOSUB  4010 S  IF  ABS(C-Cl) [ .0001  THEN  110 sC-C 1 S GOTO  90 
110  PRINTUSIHG  120 , R , H9, H9/N 

no  z  ~99.es  -eet.gge  -ft. gee 


130  NEXT  R 

140  INPUT  "DO  YOU  WANT  ANOTHER  CASE",A$Sir  AC- 
150  END 

4010  REM  *******  CD  SUBROUTINE  ******** 

4020  FOR  I-1T0  25 
4030  IF  I ]1 THEN  4050 
4040  C ( I) ”. 5 : GOTO  4070 
4050  IF  I ] 2 THEN  4070 
4060  C (I)“. 6 

4070  F(I)-2.667/C(I)*(EXP(6.6*C(I)l2)-l)-N9 
4030  IF  I [ 2 THEN  4110 

4090  C(I+1)-C(I)-F(I)*(C(I)-C(I-1))/ (F(I)-F(I 

4100  IF  AES(F(I))  [.01T1IEN  4120 

4110  NEXT  I 

4120  Cl-C(I) : RETURN 


"U"TliEH  150  s  GOTO 


D) 


3 


I 

\ 


NOMINAL  MODIFIED  REYNOLDS  NO.?  120 
VISCOSITY  CHANGE  RATIO?  .5 

NOMINAL  R(NOZ) /R(CAP)  NR'  AT  HEW  VISCOSITY  NR/NRO 


0.00 

66. 609 

0.555 

0.50 

73.509 

0.612 

1.00 

76.694 

0.  63  9 

1.50 

73.500 

0.  654 

2.00 

79.056 

0.663 

2.50 

80.456 

0.670 

3.00 

81.045 

0.075 

3.50 

SI. 495 

0.679 

4.00 

31.051 

0.  682 

4.50 

32.133 

0.034 

j 

5.00 

32.375 

0.  686 
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Term 

A 

AO 

A1 

A2 

A3 

A4 

A5 

A6 

A7 

A8 

B 

B1 

B5 

B6 

B7 

B8 

C 

CO 

C2 

C3 

C4 

C5 

0 

DO 

D1 

06 


TABLE  A-1.  COMPUTER  VARIABLE  NAME  DEFINITIONS  FOR  LAMINAR 
JET  DEVICE  PORT  CHARACTERISTIC  PROGRAM 


Definition 


Net  entrainment  coefficient/  a^ 

A0 

Term  used  in  calculation  of  Ry 

a2 

During  output  computation,  A3  is  half  the  sum  of  the  supply  and  control 
flows.  During  dynamic  computations,  A3  *  A^. 

Jet  impedance 

Total  flow  available  at  output  when  jet  is  deflected 
Switch,  1  or  0,  for  output  resistance  calculation  past  "knee" 

Output  flow  at  "knee" 

Output  pressure  at  "knee" 

Bt 

Subtrainment  radius 


sp 


sp 


cd 

Bounding  plate  loss  coefficient,  C^gp 
c9 

Control  edge  spill-back  effect  pressure  recovery  coefficient 
Limit  of  flow  for  Qc  vs  Pc  plot 

Cj 

cdbs 

Displacement  thickness 
Ay  for  Rv  integration 
Gain  limiter  -  100 
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TABLE  A-1.  COMPUTER  VARIABLE  NAME  DEFINITIONS  FOR  LAMINAR 
JET  DEVICE  PORT  CHARACTERISTIC  PROGRAM  (Cont'd) 


Term  Definition 


P4 

Phase  plot  increment 

P5 

Pr<6) 

P6 

Pr(6  -  0) 

P7 

Self-staged  output  pressure 

P8 

Phase  in  radians 

P9 

Phase  in  degrees 

C 

2e 

QO 

«o 

Q2 

AQc 

Q4 

°SB 

Q5 

Flow  recovery 

Q7 

Self-staged  output  flow 

R 

Side  vent  resistance,  Rgy 

RO 

Rc 

R1 

% 

R3 

Ri 

R4 

ROCL 

V2 

xcd 

V5 

PQ(APC  ■  -0.01)  matrix 

V6 

PQ(PC  *  0)  matrix 

V7 

PQ(PC  -  +0.01)  matrix 

W 

Bc 

HO 

£ 

c 

W3 

B°min 

H5 

0o(Pc  -  -0.01)  matrix 

W6 

Q„(PC  -  0)  matrix 

H7 

QgtPg  ■  +0.01)  matrix 
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TABLE  A-1.  COMPUTER  VARIABLE  NAME  DEFINITIONS  FOR  LAMINAR 
JET  DEVICE  PORT  CHARACTERISTIC  PROGRAM  (Cont'd) 

Term  Definition 

W9  Normalized  radian  frequency 

X  Virtual  nozzle  length 

XO  Jet  potential  core  length 

XI  Xc 

X2  Virtual  origin 

X3  XQ 

X4  Starting  value  of  Xgp 

X5  Ending  value  of  Xgp 

X6  Step  of  Xsp 

Y  Entrainment  streamline 

YO  Streamline  location 

Y1  Lateral  distance,  Y 

Y5  Running  total  of  W5 

Y6  Running  total  of  W6 


Y7 


Running  total  of  W7 
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Figure  A-l.  Printout  for  HDL  LJARS 
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Standard  silhouette  of  laminar  proportional  amplifier 
(model  3. 1.1.8) . 
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